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Abstract. We discuss the algebro-geometric initial value problem for the 
Toda hierarchy with complex-valued initial data and prove unique solvabil- 
ity globally in time for a set of initial (Dirichlet divisor) data of full mea- 
sure. To this effect we develop a new algorithm for constructing stationary 
complex- valued algebro-geometric solutions of the Toda hierarchy, which is of 
independent interest as it solves the inverse algebro-geometric spectral prob- 
lem for generally non-self-adjoint Jacobi operators, starting from a suitably 
chosen set of initial divisors of full measure. Combined with an appropriate 
first-order system of differential equations with respect to time (a substitute 
for the well-known Dubrovin equations), this yields the construction of global 
algebro-geometric solutions of the time-dependent Toda hierarchy. 

The inherent non-self-adjointness of the underlying Lax (i.e., Jacobi) oper- 
ator associated with complex-valued coefficients for the Toda hierarchy poses 
a variety of difficulties that, to the best of our knowledge, are successfully 
overcome here for the first time. Our approach is not confined to the Toda 
hierarchy but applies generally to 1 + 1-dimcnsional completely integrable (dis- 
crete and continuous) soliton equations. 



1. Introduction 

The principal aim of this paper is an explicit construction of unique global so- 
lutions of the algebro-geometric initial value problem for the Toda hierarchy with 
complex- valued initial data. More precisely, we intend to describe a solution of the 
following problem: Given p g No, assume a 1 -- ^,^ ^ to be complex- valued solutions 
of the pth stationary Toda system s-Tl p (a, b) — associated with a prescribed non- 
singular hyperelliptic curve 1C P of genus p and let r 6 No; we want to construct 
unique global solutions a = a(t r ),b = b(t r ) of the rth Tl flow Tl r (a, b) = with 
a(to,r) = a,(°\b(to,r) — for some to, r € R. Thus, we seek a unique global 
solution of the initial value problem 

Tl r (a,6) = 0, 

M)L^. p = (a< >,6<°>), ^ 
s-Tl p (a(°>,&( )) =0 (1.2) 
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for some to <r 6K, p,r€No, where a — a(n, t r ), b — b(n, t r ) satisfy 
a:Zxl^C\{0}, t:Zxl^C, 

(1.3) 

a(-,t),b(-,t) gC z , iGR, a(n, ■ ),b(n, ■) G C^R), n g Z. 

In the special case of a self-adjoint Lax (i.e., Jacobi) operator L, where a and 
6 are real-valued and bounded, the actual solution of this algebro-geometric initial 
value problem consists of the following two-step procedure discussed in detail in [3] 
(see also [H Sect. 1.3], 32, Sect. 8.3] )Q 

(i) An algorithm that constructs finite nonspecial divisors T^[i(n) G Sym p (/C p ) 
in real position for all n G Z starting from an initial Dirichlet divisor X?£( no ) G 
Sym p (/C p ) in an appropriate real position (i.e., with Dirichlet eigenvalues in appro- 
priate spectral gaps of L). "Trace formulas" of the type (|3.25p and (|3.26p then 
construct the stationary real- valued solutions a^°\ of s-TL(a, b) = 0. 

(ii) The first-order Dubrovin-type system of differential equations (|5.42[) . aug- 
mented by the initial divisor T>^ ng to \ = P/i( no ) together with the analogous "trace 
formulas" (|5.40p , (|5.4ip then yield unique global real- valued solutions a — a(t r ) , b — 
b(t r ) of the rth Tl flow Tl r (a, b) = satisfying a(t 0>r ) = a (0) , b(t , r ) = & (0) - 

This approach works perfectly in the special self-adjoint case where the Dirichlet 
divisors fi(n,t r ) = (/ti(n, t r ), . . . ,fx p (n,t r )) G Sym p (AC p ), (n,t r ) G Z x R, yield 
Dirichlet eigenvalues /ii(n, t r ), . . . , fJ- p (n, t r ) of the Lax operator L situated in p 
different spectral gaps of L on the real axis. In particular, for fixed (n, t r ) G Z x R, 
the Dirichlet eigenvalues fij(n,t r ), j = 1, . . . ,p, are pairwise distinct and formulas 
(|5.4ip for a and (|5.42p for (d/dt r )fj.j(n,t r ), j — 1, . . . ,p, are well-defined. 

This situation drastically changes if complex- valued initial data a^°\b^ ' or 
p,(n ,t r ) are permitted. In this case the Dirichlet eigenvalues /j,j(n,t r ), j — 
1, . . . ,p, are no longer confined to well separated spectral gaps of L on the real axis 
and, in particular, they are in general no longer pairwise distinct and "collisions" 
between them can occur at certain values of (n,t r ) G Z x R. Thus, the stationary 
algorithm in step (i) as well as the Dubrovin-type first-order system of differential 
equations (|5.42p in step (ii) above, breaks down at such values of (n, t r ). A priori, 
one has no control over such collisions, especially, it is not possible to identify ini- 
tial conditions T^p,(n 0> t r) a ^ sorne (^Ojio.r) G Z x R which avoid collisions for all 
(n, t) G Z x R. We solve this problem head on by explicitly permitting collisions in 
the stationary as well as time-dependent context from the outset. In the stationary 
context, we properly modify the algorithm described above in step (i) in the self- 
adjoint case by alluding to a more general interpolation formalism (cf. Appendix |B|) 
for polynomials, going beyond the usual Lagrange interpolation formulas. In the 
time-dependent context we replace the first-order system of Dubrovin-type equa- 
tions (|5.42p . augmented with the initial divisor P/}( no> t r ), by a different first-order 
system of differential equations (|6.27p with initial conditions (|6.28p which focuses 
on symmetric functions of fj,i(n,t r ), . . . , fi p (n,t r ) rather than individual Dirichlet 
eigenvalues fij(n,t r ), j = 1, . . . ,p. In this manner it will be shown that collisions 
of Dirichlet eigenvalues no longer pose a problem. 

In addition, there is a second nontrivial complication in the non-self- adjoint 
case: Since the Dirichlet eigenvalues fj,j(n,t r ), j = l,...,p, are no longer confined 
to spectral gaps of L on the real axis as (n,t r ) vary in Z x R, it can no longer 



We freely use the notation of divisors of degree p as introduced in Appendix fA] 
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be guaranteed that fij(n,t r ), j = 1, . . . ,p, stay finite for all (n,t r ) EZxI. As 
discussed in Section [4] in the stationary case, this phenomenon is related to certain 
deformations of the algebraic curve IC P under which for some uq £ Z, a(rio) — > 
and Hj(riQ + 1) — > oo for some j G {1, . . . ,p}. We solve this particular problem 
in the stationary as well as time-dependent case by properly restricting the initial 
Dirichlet divisors D^ n(uto r j S Sym p (/C p ) to a dense set of full measure. 

Summing up, we offer a new algorithm to solve the inverse algebro-geometric 
spectral problem for generally non-self-adjoint Jacobi operators, starting from a 
properly chosen dense set of initial divisors of full measure. Combined with an 
appropriate first-order system of differential equations with respect to time (a sub- 
stitute for the well-known Dubrovin equations), this yields the construction of global 
algebro-geometric solutions of the time-dependent Toda hierarchy. 

We emphasize that the approach described in this paper is not limited to the 
Toda hierarchy but applies universally to constructing algebro-geometric solutions 
of 1+1 -dimensional integrable soliton equations. In particular, it applies to differential- 
difference (i.e., lattice) systems and we are now in the process of applying it to the 
Ablowitz-Ladik hierarchy. Moreover, the principal idea of replacing Dubrovin- type 
equations by a first-order system of the type (|6.27p is also relevant in the context of 
general non-self-adjoint Lax operators for the continuous models in 1 + 1-dimensions. 
(In particular, the models studied in detail in [13j can be revisited from this point 
of view.) We also note that while the periodic case with complex-valued a, b is of 
course included in our analysis, we throughout consider the more general algebro- 
geometric case (in which a,b need not even be quasi-periodic). 

Finally we briefly describe the content of each section. Section [2] presents a quick 
summary of the basics of the Toda hierarchy, its recursive construction, Lax pairs, 
and zero-curvature equations. The stationary algebro-geometric Toda hierarchy so- 
lutions, the underlying hyperelliptic curve, trace formulas, etc., are the subject of 
Section [3j A new algorithm solving the algebro-geometric inverse spectral problem 
for generally non-self-adjoint Jacobi operators is presented in Section^] In Section^] 
we briefly summarize the properties of algebro-geometric time-dependent solutions 
of the Toda hierarchy and formulate the algebro-geometric initial value problem. 
Uniqueness and existence of global solutions of the algebro-geometric initial value 
problem as well as their explicit construction are then presented in our final and 
principal Section [6j Appendix [X] reviews the basics of hyperelliptic Riemann sur- 
faces of the Toda-type and sets the stage of much of the notation used in this 
paper. Various interpolation formulas of fundamental importance to our station- 
ary inverse spectral algorithm developed in Section U are summarized in Appendix 
iBl Finally, Appendix [C] summarizes asymptotic spectral parameter expansions of 
various quantities fundamental to the polynomial recursion formalism presented in 
Section [21 These appendices support our intention to make this paper reasonably 
self-contained. 

2. The Toda Hierarchy in a Nutshell 

In this section we briefly review the recursive construction of the Toda hierarchy 
and associated Lax pairs and zero-curvature equations following [BJ, [141 Sect. 1.2], 
and [51 Ch. 12]. 

Throughout this section we make the following assumption: 
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Hypothesis 2.1. Suppose 

a, b £ C Z and a{n) ^ for all n £ Z. (2.1) 

Here C J denotes the set of complex-valued sequences indexed by J C Z. 
We consider the second-order Jacobi difference expression 

L = aS+ +a-S- +b, (2.2) 

where S denote the shift operators 

(S ± f)(n) = /±(n) = /(nil), n £ Z, f £ C z . (2.3) 

To construct the stationary Toda hierarchy we need a second difference expres- 
sion of order 2p + 2, p £ No, defined recursively in the following. We take the 
quickest route to the construction of P2 P +2, and hence to the Toda hierarchy, by 
starting from the recursion relations (|2.4p - (|2.6p below. 

Define {fe}i<EN and {gi}i e n recursively by 

/o = l, g = - Cl , (2.4) 

2//+1+& + gj -2bft = 0, ££N , (2.5) 

g e +i-gi +1 + 2(a 2 f+ -(a-) 2 fi)-b(g e -gi)=0, I £ N . (2.6) 
Explicitly, one finds 



(2.7) 



= -2a 2 (6 + + 6)+ci(-2a 2 )-c 3 , etc. 

Here {c^jigN denote undetermined summation constants which naturally arise when 
solving (lf4l- (l2~B)) . 

Subsequently, it will also be useful to work with the corresponding homogeneous 
coefficients fj and gj, defined by vanishing of the constants Cfc, k £ N, 



fo = 


1, 


fl = 


6 + ci, 


h = 


a 2 + (eT 


go = 


-ci, 


9i = 


-2a 2 - c 


92 = 


-2a 2 (6+ 



(2.8) 



Hence. 



introducing 



/e = ct-kfk, gt = ^ Cf - fc 5fe ~ c «+i' ^ G N o, (2.9) 



fe=0 fe=l 



c = 1. (2.10) 
Next we define difference expressions i"2p+2 of order 2p + 2 by 

P 2p+2 = -L p+1 + ^( 5 ^ + 2 a /^+)L^ +f p+1 , P eN . (2.11) 
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Introducing the corresponding homogeneous difference expressions P 2p +2 defined 

by 

^+3 = ^+2^=0,^,...,,, ^N„, (2.12) 

one finds 

p 

P2p+2 = C p -eP2l+2- (2-13) 
1=0 

Using the recursion relations (I2.4[) — (|2.6[) . the commutator of P2 P +2 and L can be 
explicitly computed and one obtains 

[P 2p+2 , L] = - a( 9 ; + g p + f+ +1 + f p+1 - 2b+f+)S+ 

+ 2( - b(g p + f p+1 ) + a 2 f+ - (a-) 2 f~ + b 2 f p ) 

- ^ (g P + g p + f P +i + f- +1 - 2bf p )S~, p E N . (2.14) 

In particular, (L,p2 P +2) represents the celebrated Lax pair of the Toda hierarchy. 
Varying p e No, the stationary Toda hierarchy is then defined in terms of the 
vanishing of the commutator of P2 P +2 and L in (|2.14p . that is, 

[P 2p+2 ,L] =s-Tl p (a,6) -0, peN . (2.15) 

Thus one finds 

9p + 9p + fp+i + fp+i - 2bf P = 0, (2.16) 
-K9p + fp+i) + a 2 f+ - (a-) 2 f- + b 2 f p = 0. (2.17) 
Using (|2.5|) with j = p one concludes that 12 . 16[) reduces to 

fp+i = /p+i, (2-18) 

that is, fp+i is a lattice constant. Similarly, one infers by subtracting b times (|2.16p 
from twice (|2.17j) and using (|2.6p with j — p, that g p +i is a lattice constant as well, 
that is, 

9p+i=9 p +i- (2-19) 
Equations (|2.18|) and (|2.19p give rise to the stationary Toda hierarchy, which is 
introduced as follows 

s-Tl p (a, b) = ( f P +1 ~ fp + 1 ) =0, pe N . (2.20) 
\9p+i-9p+iJ 

Explicitly, 

b+ -b 

,2((c 



s-Tl (a,6) 



sTU(ab)-( (a + ) 2 -(a-f + (b + ) 2 -b 2 \ 

s-ll^a, b) - ^ 2(a - )2(6 + b - ) _ 2a 2 {b + +b) J (2-21) 

( b+ - b \ n 

+ Cl {2((a-) 2 -a 2 )) =°' 6tC -' 

represent the first few equations of the stationary Toda hierarchy. By definition, 
the set of solutions of (|2.20p . with p ranging in No and c< 6 C, f € N, represents 
the class of algebro-geometric Toda solutions. 

In the following we will frequently assume that a, b satisfy the pth stationary 
Toda system. By this we mean it satisfies one of the pth stationary Toda equations 
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after a particular choice of summation constants ct 6 C, I = 1, . . . ,p, p £ N, has 
been made. 

In accordance with our notation introduced in (12. 8|) and (|2.12|) , the corresponding 
homogeneous stationary Toda equations are defined by 

s-fi p (a,6)=s-Tl p (a,6)| cf=0;fclj p = 0, p e N„. (2.22) 

Next, we introduce polynomials F p and G p +\ of degree p and p+ 1, with respect 
to the spectral parameter z 6 C by 

p p 

F p (z) = J2 fp-^ = c p-^( z )> ( 2 - 23 ) 

1=0 £=0 

p p+1 

G p+1 (z) = -z p+1 + Y,9 P -iz e + f p+1 = Y,Cp+i-zGe{z) (2.24) 

with Fi and G^ denoting the corresponding homogeneous polynomials defined by 
F (z) = F Q (z) = 1, 



Fi(z) = F t (z)\ Ck=Qtk=lt ...,, = ^ /^ fe , ^N„, (2.25) 

fc=0 

G (z) = G (z)| ci=0 = 0, Gi(z) = Gi(z) = -z - 6, 
G, +1 (z) - G, +1 (z)| Cfe=Q fc=1 ^ = + ]T^_ fe z fe + /,+!, t G N. (2.26) 

A;=0 

Explicitly, one obtains 
F = l, 

Fx = z + b + ci, 

F 2 = z 2 + bz + a 2 + {a-f + b 2 + ci(z + 6) + c 2 , etc., (2.27) 
Gx = -z + b, 

G 2 = -z 2 + (a - ) 2 - a 2 + b 2 + a(-z + b), etc. 

Next, we study the restriction of the difference expression P 2p +2 to the two- 
dimensional kernel (i.e., the formal null space in an algebraic sense as opposed to 
the functional analytic one) of (L — z). More precisely, let 

ker(L -z) = {^:Z-»CU {oo} | (L - z)tp = 0}. (2.28) 

Then (f2TTT]) implies 

P2 P+ 2 Wer(L-z)= {2aF p (z)S+ + G p+1 {z))\^ (L _ z) . (2.29) 

Therefore, the Lax relation (|2.15p becomes 

2(z - b+)F+ - 2(z - b)F p + G+ +1 - G- +1 = 0, (2.30) 

2a 2 F+ - 2{a-) 2 F p + {z - b){G p+1 - G p+1 ) = 0. (2.31) 

Additional manipulations yield 

2( 2 - b)F p + G p+1 + G p+l = 0, (2.32) 

{z - b) 2 F p + {z- b)G p+1 + a 2 F+ - {a~) 2 F p = 0. (2.33) 



THE TODA HIERARCHY INITIAL VALUE PROBLEM 



7 



Indeed, adding G p +i — G p +i to the left-hand side of (12.30|) (neglecting a trivial 
summation constant) yields (|2.32|) and inserting (|2.30[) into (|2.31|1 then implies 
(|2.33p . Varying p £ No , equations (12.32|) , (|2.33[) provide an alternative description 
of the stationary Toda hierarchy. 

Combining equations (|2.31|) and (|2.32[) one concludes that the quantity 

R2 V +2{z) = G p+1 {z,n) 2 -4a(n) 2 F p (z,n)F+(z,n) (2.34) 

is a lattice constant, and hence depends on z only. Thus, one can write 

2J9+1 

R2 P+ 2(z) = J] (z - E m ), {E m } 2 n p +l C C. (2.35) 

m— 

One can decouple (I2.32p and (]2 .33[) to obtain separate equations for F p and G p +\. 
For instance, computing G p +i from (|2 . 33[) and inserting the result into (|2.32|) yields 
the following linear difference equation for F p 

(z - b) 2 (z - b~)F p -(z- b-) 2 (z - b)F~ + ((a-fF- - a 2 F+)(z - b~) 

+ {(a—) 2 F p - - (a-) 2 F p )(z -b) = 0. (2.36) 

Similarly, insertion of (|2.33|) into (|2.34|) permits one to eliminate G p +i and results 
in the following nonlinear difference equation for F p , 

(z - bfF 2 - 2a 2 (z - b) 2 F p F+ - 2( fl -) 2 (z - b) 2 F p F~ + a 4 (F+) 2 

+ (a-)\F-) 2 - 2a 2 {a-) 2 F+F p = (z - b) 2 R 2p+2 (z). (2.37) 

On the other hand, computing F p in terms of G p +\ and G p+1 using (|2.32|) and 
inserting the result into (I2.33|) yields the following linear difference equation for 
G P +i 

a 2 (z - b~)(G+ +1 + G p+1 ) - (a-) 2 (z - b+)(G p+1 + G^) 

+ (z- b-)(z -b)(z- b+)(G- +1 - G p+1 ) = 0. (2.38) 

Finally, inserting the result for F p into (|2.34[) yields the following nonlinear differ- 
ence equation for G p +i 

(z - b)(z - b+)G 2 p+1 - a 2 (G p+1 + G P+1 )(G P+1 + G+ +1 ) 

= (z-b)(z-b+)R 2p+2 . (2.39) 

Equations (|2.37[) and (|2.39|) can be used to derive nonlinear recursion relations for 
the homogeneous coefficients ft and gg (i.e., the ones satisfying l|2.8[) in the case of 
vanishing summation constants) as proved in Theorem IC.ll in Appendix [Cj This 
has interesting applications to the asymptotic expansion of the Green's function of 
L with respect to the spectral parameter. In addition, as proven in Theorem lC.il 
(|2.37[) leads to an explicit determination of the summation constants c\ , . . . , c p in 

s-Tl p (a,6)=0, peN , (2.40) 

in terms of the zeros Eq, . . . , E 2p +i of the associated polynomial i?2p+2 in (|2.35[) . 
In fact, one can prove (cf. Theorem I C.l[) that 

c k =c k {E), k = l,...,p, (2.41) 
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where 

k 

„ / p \ _ _ ( 2 io) ! • ■ ■ ( 2 j2p+l)! pj pJ2p+ i 

~ ,o,..i +1 =o 22 "^ ! ) 2 • ' ' ^v + mih -!)••• (2j 2p+1 - l)^ ' " ' 

30 H hj2p+l=fc 

fc = l,...,p, (2.42) 

are symmetric functions of E = (Eg, . . . , E 2p+ i). 

We emphasize that the result (|2.29[) is valid independently of whether or not 
P 2p + 2 and L commute. However, the fact that the two difference expressions P 2p +2 
and L commute implies the existence of an algebraic relationship between them. 
This gives rise to the BurchnalL-Chaundy polynomial for the Toda hierarchy first 
discussed in the discrete context by Naiman [28], [29] , 

Theorem 2.2. Assume Hypothesis 12.11 fix p € No and suppose that P 2p +2 and L 
commute, [P 2p+2l L] = 0, or equivalently, assume that s-Tl p (a, b) = 0. Then L and 
P2 P +2 satisfy an algebraic relationship of the type (cf. ()2.35j) ) 

Fp(L, P2 P +2) = Pip+2 ~ R2 P +2{L) = 0, 

2 p± x (2.43) 
R2 P +2(z)= Y[{z-E m ), zeC. 

The expression J- p (L, P2 P +2) is called the Burclmall-Chaundy polynomial of the 
Lax pair (L, P 2p+2 ) and it will be used in Section [3] to introduce the underlying 
hyperelliptic curve associated with the stationary Toda system s-Tl p (a, b) — (cf. 

©3D)- 

Next we turn to the time-dependent Toda hierarchy. For that purpose the func- 
tions a and b are now considered as functions of both the lattice point and time. 
For each equation in the hierarchy, that is, for each p, we introduce a deforma- 
tion (time) parameter t p g R in a, b, replacing a(n),b(n) by a(n,t p ),b(n,t p ). The 
second-order difference expression L (cf. (|2.2p ) now reads 

L(t p ) = a( • , t p )S + + a~( ■ , t p )S- + b{ • , t p ). (2.44) 

The quantities {ft}i^n , {ge}eeN , and P2 P +2, p6N are still defined by (|2.4|l — (|2.6[) 
and (|2.11[) . respectively. The time-dependent Toda hierarchy is then obtained by 
imposing the Lax commutator equations 

L tp (t p )-[P 2p+2 (t P ),L(t p )] = 0, t p G R, (2.45) 

varying p € No . Relation (12.45)) implies 

(a tp + a(g+ + g p + f+ +1 + f p+1 - 26+/+)) S+ 

- (-b tp + 2( - b(g p + f p+1 ) + a 2 f+ - (a-) 2 f- + b 2 f p )) (2.46) 

+ (a tp + a(g+ + g p + f+ +1 + f p+1 - 26+/+))" S~ = 0. 
Applying the same method we used to derive (|2.18[) and (|2.19[) one concludes 
= L tp -[P 2p+2 ,L] 
= (a tp - a(/+ +1 - / p+ i)) 5+ - (-b tp - gp+i + g p+1 ) 

+ (at p - a(f+ +1 - f p+1 )y S- . (2.47) 
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Varying p <E No , the collection of evolution equations 

TL (a, b) = ( a \ ~ a(/ P ++1 ~ /p+l) ^) = 0, (n, *-) G Z x R, p e No (2.48) 
V % + 9 P +i ~ ffp+i / 



then defines the time-dependent Toda hierarchy. Explicitly. 



Tl (n h\-( ata ~ a ( 6+ ~ 6 ) 



= 0, 



Tl (a b) N ~ <^ - (a ^ 2 + ^ b ')\ (2 49) 

Tlx (a, b) - ^ ^ + 2 v (a _ )2(6 + b _ } _ 2a2(b+ + h >j (2.49) 

( -a{b+-b) \ n 
+ Cl v -2(a 2 -(a-) 2 )J =0 ' eta ' 

represent the first few equations of the time-dependent Toda hierarchy. The system 
of equations, Tlo(a, b) = 0, is of course the Toda system. 

The corresponding homogeneous Toda equations obtained by taking all summa- 
tion constants equal to zero, ci = 0, I = 1, . . . ,p, are then denoted by 

T\ p (a,b) = T\ p (a,b)\ ct=oe=h p . (2.50) 

Restricting the Lax relation (|2.45[) to the kernel ker(L — z) one finds that 

= (L tp [P 2p+2 , L])\ kcr{L _ z) = (L tp + (L - z)P 2p+2 ) | kcr(i _ z) (2.51) 

a(^-^ + 2(z- b+)F+ - 2{z - b)F p + G+ +1 - G~ +1 ) S+ 



+ (b tp + (z~ b)°^ + 2(a-) 2 F p - 2a 2 F+ 
+ ( z -b)(G- +1 -G p+1 ))\ 



(2.52) 



kcr(L- 



a 


a~ 


= -2(z 


K = 


-(* 


cu 
- b)-^- 

a 



Hence one obtains 

^ = -2(z - b+)F+ + 2{z - b)F p + G p+l - G+ +1 , (2.53) 

7 + 2a 2 F+ -(z- b)(G p+1 - Gp +1 ). (2.54) 

Further manipulations then yield, 

a tp = -a(2(z - b+)F+ + G+ +1 + G p+1 ), (2.55) 

b tp =2((z-b) 2 F p + {z-b)G p+1 +a 2 F+ -(a~) 2 F p ). (2.56) 

Indeed, (|2.55|) follows by adding G p +i — G p +i to (|2.53|) (neglecting a trivial summa- 
tion constant), and an insertion of (|2.55[) into (|2.54p implies (|2.56l) . Varying p G No, 
equations (|2.55p and (|2.56[) provide an alternative description of the time-dependent 
Toda hierarchy. 

Remark 2.3. From (|2~4|) - ([2Hj) and ([2^23]) . (|2~24|) one concludes that the coeffi- 
cient a enters quadratically in F p and G p+ i, and hence the Toda hierarchy (I2.48[) 
(respectively (|2.20[1 ) is invariant under the substitution 

a -> a £ = {e(n)a(n)} ne z, e{n) G {1, -1}, n G Z. (2-57) 
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We conclude this section by pointing out an alternative construction of the Toda 
hierarchy using a zero-curvature approach instead of Lax pairs (L, p2 P +2)- To this 
end one defines the 2x2 matrices 



U(z) = 



1 

-aT I a (z — b) I a 



V,, (z) - \ G p+^ 
Vp+l[Z) ~ ^-2a-F p (; 



2a~F~(z) 
z) 2(z-b)F p {z) + G p+1 (z) 



P£N„. 



(2.58) 
(2.59) 



Then the stationary part of this section can equivalently be based on the zero- 
curvature equation 







uv p+1 - V+ +1 U 

I 

r((z-b+)F+ -(z_-b)F p 
+ 2 1 (C p +i - Gp+i)) 



(2.60) 







V 



a 2 F+-( a -) 2 F p 
+2- 1 (z-6)(G p+1 -G+ +1 
+ {z-b) 2 F p -(z-b+){z- b)F+J 



Thus, one obtains ([2~30]) from the (2, l)-entry in ([2760]) . Insertion of (|2~30| into the 
(2,2)-entry of ((27601) tnen yields ([275T]) . Thus, one also obtains (|2~32]) and hence 
the (2, 2)-entry of V p +\ in (|2.59p simplihes to 



i(z) = -G +1 (z) 



(2.61) 



in the stationary case. Since det(C/(z,n)) = a~(n)/a(n) / 0, n £ Z, the zero- 
curvature equation (|2 . 60[) yields that det(V p +i (z,n)) is a lattice constant (i.e., in- 
dependent of n G Z). The Burchnall-Chaundy polynomial F p (y,z) (cf. (|2.43|) and 
especially, the hyperelliptic curve (|3.1j) ) is then obtained from the characteristic 
equation of V p +i{z) by 

det(y/ 2 - Vp+i(z,n)) 
= y 2 + dct(Vp + i(z,n)) 

= y 2 - Gpl^z,™) 2 + 4oT (n) 2 F p (z 7 n)F p {z, n) 
' R 2p+ 2{z) = 0, 



y 



(2.62) 



using (|2.61|) . (Here I2 denotes the identity matrix in C 2 .) Similarly, the time- 
dependent part (|2.44p - (|2.56p can equivalently be developed from the zero-curvature 
equation 



Q = U t + UK 



P+i 



I 



V 



(2.63) 
\ 



a (( a t P / a ) - ( a t p / a )) 
-a- (2(z - b+)F+ - 2{z - b)F p 

+ (G p+ i - G p+1 )) 





-b tp - {z-b){a t Ja) 
+2a 2 F+-2(a-) 2 F p 
+(z-b)(G v+1 -G+ +1 ) 
+2{z - b) 2 F p - 2{z - b+){z - b)F+ j 



The (2, l)-entry in ([2~53"j) yields ([2331) . and inserting (|2~53"1) into the (2, 2)-entry of 
(I2.63|) yields (|2.54p and hence also the basic equations defining the time-dependent 
Toda hierarchy in (j2"35)) . (j2~56l) . 
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3. Properties of Stationary Algebro-Geometric Solutions 
of the toda hierarchy 

In this section we present a quick review of properties of algebro-geometric so- 
lutions of the stationary Toda hierarchy. Since this material is standard we omit 
all proofs and just refer to [6] (cf. also [HI Sect. 1.3], [311 Chs. 8, 9]) for detailed 
presentations and an extensive list of references to the literature. 

For the notation employed in connection with elementary concepts in algebraic 
geometry (more precisely, the theory of compact Riemann surfaces), we refer to 
Appendix [AJ 

Returning to Theorem l2.21 we note that (|2.43[) naturally leads to the hyperelliptic 
curve IC P of genus p 6 No, where 

Kp : T p (z, y) = y 2 - R 2p +2(z) = 0, 

2p+l ]\ 

R 2p+2 (z)=l[(z-E m ), {EJScC. 

m— 

Throughout this section we make the following assumption: 

Hypothesis 3.1. Suppose that 

a, b e C z and a(n) ^ for all n G Z. (3.2) 

In addition, assume that the hyperelliptic curve JC P in (|3.ip is nonsingular, that is, 
suppose that 

E m ^ E m i for m 7^ m! , m, m = 0, . . . , 2p + 1. (3-3) 

The curve K, p is compactified by joining two points Poo± , Poc + 7^ Poo- , at infinity. 
For notational simplicity, the resulting curve is still denoted by JC p . Points P on 
fcp\{Poo + 1 Poo- } are represented as pairs P = (z, y), where y( ■ ) is the meromorphic 
function on K, p satisfying T p (z, y) — 0. The complex structure on IC P is then defined 
in the usual way, see Appendix [Al Hence, K, p becomes a two-sheeted hyperelliptic 
Riemann surface of genus p 6 No in a standard manner. 

We also emphasize that by fixing the curve K, p (i.e., by fixing Eq, ■ ■ ■ ,i?2p+i), 
the summation constants c\ , . . . , c p in the corresponding stationary s-Tl p equation 
are uniquely determined as is clear from (|2.41[) and (|2.42[) . which establish the 
summation constants cj, as symmetric functions of Eq, . . . , i?2 P +i- 

For notational simplicity we will usually tacitly assume that pel. The trivial 
case p = 0, which leads to a(n) 2 = {E x - S ) 2 /16, b(n) = (Eq + E{)/2, n E Z, is of 
no interest to us in this paper. 

In the following, the zero:Q of the polynomial F p (- ,n) (cf. (|2.23p ) will play a 
special role. We denote them by {^j{n)} p - =1 and write 

p 

F p {z,n) = Y[{z- l i J {n)). (3.4) 

j=l 

The next step is crucial; it permits us to "lift" the zeros Hj of F p from C to the 
curve /C p . From (|2.34ll and (|3.4|) one infers 

R2 P +2(z) - G p+ i(z) 2 =0, z <E {^j,^t}j,k=i,..., P - (3.5) 

2 If a, b S £°°(Z), these zeros are the Dirichlet eigenvalues of a bounded operator on £ 2 (Ij) 
associated with the difference expression L and a Dirichlet boundary condition at n £ Z. 



12 



F. GESZTESY, H. HOLDEN, AND G. TESCHL 



We now introduce {fij(n)}j=i,...,p C IC P by 

p,j(n) = (jMj(ri),-Gp +1 (fj,j(n),ri)) E 1C P , j = 1, ...,p, n E Z. (3.6) 

Next, we recall equation (12.34| and define the fundamental meromorphic function 
<f>(-,n) on IC P by 

mn) = y- G p+f^) (3 .7) 

-2a(n)F p (z,n+ 1) 

P = (z,y) e K p , n e Z, 
with divisor (0( • , n)) of </>(•, n) given by 

(<H ',»))= ^Poo+^n+l) ~ T> P^_Kn), (3-9) 

using (|3.4p and (|3.6p . Here we abbreviated 

A = {Ai: ■ ■ ■ j A P } G Sym p (/C p ) (3.10) 

(cf. the notation introduced in Appendix [XJ . We note that several fij (n) may be 
equal for a given lattice point neZ. Moreover, since — G p +i(jUj(n), n) takes on the 
same value for all coinciding zeros Hj(n), no finite special divisors ~Dji(n) can ever 
arise in (cf. also Lemma I3T4")) . 

The stationary Baker-Akhiezer function , n, no) on K, p \ {Poo± } is then de- 
fined in terms of <f>( ■ , n) by 

frciLo^"*) for »>«o+ii 

4>(P, n, no) = < 1 for n = no, (3-11) 

llKi 1 ^.™)" 1 f or n<n -l, 
PG/CpXfP^}, (n,n ) GZ 2 , 

with divisor [ip( ■ , n, no)) of ip(P, n, no) given by 

(ip( ■ , n, n )) = P^ ( „) - Z>A(„ ) + (n - n )(2?p =o+ - X>p^ ). (3.12) 

For future purposes we also introduce the following Baker-Akhiezer vector, 

*(P,n, no) = (^^f) ' P G ^ {Foo±} ' ( "' no) G (3 - 13) 

Basic properties of <J), ip, and ^ are summarized in the following result. We 
abbreviate by 

W(f,g) = a(fg + -f+g) (3.14) 

the Wronskian of two complex-valued sequences / and g, and denote P* — (z, —y) 
for P = (z,y) G tCp. 

Lemma 3.2. Assume Hupothesis YS.W and suppose that a, b satisfy the pth stationary 
Toda system (|2.20p . Moreover, let P = (z,y) E IC p \{Poo ± } and (n,no) E Z 2 . Then 
<j) satisfies the Riccati-type equation 

a<i>{P) + a~ <p- {Py 1 = z-b, (3.15) 
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(3.22) 



as well as 

cj)(P)<j)(P*) = 2£1 (3.16) 

« p ) + « p ') = -^ <" 7 > 

* (P) -« P *> = SW (3 - 18) 

Moreover, ip and VP satisfy 

(L - z(P))ip(P) = 0, (P 2p+2 -y(P))ip{P) = 0, (3.19) 

*+(P) = U(zMP), yV(P)=V p+1 V(P), (3.20) 

V(P, n, n )V(P*, n, no) = =^ ' ^ (3.21) 

Fp(z,n ) 

a(n)(ip(P,n,n )ip(P*,n + l,n ) + Tp(P* ,n,n )ip(P,n + l,n )) 
_ G p+ i(z,n) 
F p (z,n ) 

W(tP(P, • , n ), • , n )) = ~ V ( 3 ' 23 ) 

P P (z, n ) 

Combining the polynomial recursion approach with (|3.4|) readily yields trace 
formulas for the Toda invariants, which are expressions of a and b in terms of the 
zeros Hj of F p . We introduce the abbreviation, 

2p+l p 
m— j—1 

Lemma 3.3. Assume Hvvothesis Yd.W and suppose that a, b satisfy thepth stationary 
Toda system (|2.20[) . Then, 

2p+l p 
m— 7=1 

in addition, if for all n £ Z, f^j(n) ^ /or j ^ k, j,k = 1, . . . ,p } then, 

«H 2 = 5 E II ( n ) - MfeH)- 1 + J (6 (2) (n) - 6(n) 2 ) , «£l (3.26) 

The case where some of the fij coincide in (|3.26|) requires a more elaborate 
argument that will be presented in Section 3J 

Since nonspecial Dirichlet divisors T)^ and the linearization property of the Abel 
map when applied to will play a fundamental role later on, we also recall the 
following facts. 

Lemma 3.4. Assume Hvvothesis Yd .11 and suppose that a, b satisfy the pth stationary 
Toda system (|2.20|) . Let T>f L , fi — {/ti, . . . , fj, p } £ Sym p (/C p ), be the Dirichlet divisor 
of degree p associated with a, b defined according to (|3.6p . that is, 

Ai( n ) = (Mi( n )i -Gp+itftW^)) G ^P) .7 = 1, • ■ ■ n S Z. (3.27) 
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Then Dp,( n ) * s nonspecial for all n e Z. Moreover, the Abel map linearizes the 
auxiliary divisor in the sense that 

fiQo(%«)) = «Qo( P A(«o)) ~ i n ~ n o)Ap ao _ (P°°+), (3-28) 
where Qo £ IC P is a given base point. 

If in addition, a, 6 € £°°(Z), then there exists a constant > such that 

M«)|<C M , j = l,...,p, neZ. (3.29) 

Remark 3.5. We note that by construction, the divisors D^^a, n € Z, as intro- 
duced in (|3.6|) are a^Z finite and hence nonspecial by Lemma \3.4\ On the other hand, 
as we will see in the next Section® given a nonspecial divisor Z?^(„ ) ; the solution 
Dfafn) of eguation (13.28[) may cease to be a finite divisor at some n£Z. 

4. An Algorithm for Solving the Inverse Algebro-Geometric 
Spectral Problem for (Non-self-adjoint) Jacobi Operators 

The aim of this section is to derive an algorithm that enables one to construct 
algebro-geometric solutions for the stationary Toda hierarchy for complex-valued 
initial data. Equivalently, we offer a solution of the inverse algebro-geometric spec- 
tral problem for general (non-self-adjoint) Jacobi operators, starting with initial 
divisors in general complex position. 

Up to the end of Section [3] the material was standard (see [S] and P31 Sect. 1.3], 
[521 Chs. 8, 9] for details) and based on the assumption that a,i £ C z satisfy the 
pth stationary Toda system (|2.20p . Now we embark on the corresponding inverse 
problem consisting of constructing a solution of (I2.20[) given certain initial data. 
More precisely, we seek to construct solutions a, b e C z satisfying the pth stationary 
Toda system (|2.20j) starting from a properly restricted set Ado of finite nonspecial 
Dirichlet divisor initial data D^r no \ at some fixed no € Z, 

pt(n ) = {£i(n ), ■ • • , Ap(™o)} G M , M c Sym p (/C p ), 
— (4.1) 
Aifao) = w( n o),-Gp+i(ti>j(no),n )), j = l,...,p. 

Of course we would like to ensure that the sequences obtained via our algorithm 
do not blow up. To investigate when this happens, we study the image of our 
divisors under the Abel map. The key ingredient in our analysis will be (|3.28[) 
which yields a linear discrete dynamical system on the Jacobi variety J(/C p ). In 
particular, we will be led to investigate solutions "Dp, of the discrete initial value 
problem 

aQo^Ato) = «Q ( p A(«o)) - (n - n )A Pao (P x+ ), 

(4.2) 

A(«o) = {Aifao), ■ . . ,/t p (n )} G Sym p (/C p ), 

where Qo G IC P is a given base point. Eventually, we will be interested in solutions 
Vp of (|4.2p with initial data 2?/t(n ) satisfying (|4.1|) and A4o to be specified as in 
(the proof of) Lemma 14.21 

Before proceeding to develop the stationary Toda algorithm, we briefly analyze 
the dynamics of (|4.2[) . 

Lemma 4.1. Let Dp,(n) be defined via (14. 2|) for some divisor T>^ no ^ 6 Sym p (/C p ). 
(i) If D^n) * s finite and nonspecial and T>^ n+1 ^ is infinite, then £>^( n +i) contains 
Poo+ bvl not P OQ _ . 
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(ii) IfDjx( n ) is nonspecial and 2?/i( n +i) is special, then 'Dp.in) contains P 00+ at least 
twice. 

Items (i) and (ii) hold if n + 1 is replaced by n — 1 and Poa + by Poo- ■ 

Proof, (i) Suppose one point in T>pr n +x) equals Poo„ and denote the remaining ones 
by £>£( n+1) . Then gJJ tells us a Po (V^ n+1 -)) + A Po (P 00+ ) = ap (£>#(„)). Since we 
assumed £>/}( ra ) to be nonspecial, we have T>fi( n ) = ^( n +i) + ^- > p o + contradicting 
finiteness of V^ n y 

(ii) We choose Pq to be a branch point such that A P (P*) = —A P (P). In par- 
ticular, if D^n+i) is special, then it contains a pair of points (Q,Q*) whose con- 
tribution will cancel under the Abel map, that is, QpJI'^n+i)) = apj^jn+i)) 
for some V^n+i) £ Sym p ~ 2 (/C p ). But invoking (|4.2p shows that ap (£>^( n )) = 
a Po (2?£(„ +1 )) + 2A Pa (P 00+ ). As Dp,( n ) was assumed nonspecial, this shows that 
£>A(n) = %(™+i) + 2X'p oc+ , as claimed. □ 

This yields the following behavior of T)[i(n) if we start with some nonspecial 
finite initial divisor V^^y. As n increases, D^n) stays nonspecial as long as it 
remains finite. If it becomes infinite, then it is still nonspecial and contains Poo + 
at least once (but not Poo-)- Further increasing n, all instances of Poo + will be 
rendered into Poo_ step by step, until we have again a nonspecial divisor that has 
the same number of Poo- as the first infinite one had Poo + - Generically, we expect 
the subsequent divisor to be finite and nonspecial again. 

Next we show that most initial divisors are nice in the sense that their iterates 
stay away from infinity. Since we want to show that this set is of full measure, it 
will be convenient for us to identify Sym p (/C p ) with the Jacobi variety J(IC P ) via 
the Abel map and take the Haar measure on J(/C p ). Of course, the Abel map is 
only injective when restricted to the set of nonspecial divisors, but these are the 
only ones we are interested in. 

Lemma 4.2. The set A4q C Sym p (/C p ) of initial divisors T)^ nf ^ for which T)^ n y 
defined via (|4.2p . is finite and hence nonspecial for all n eZ, forms a dense set of 
full measure in the set Sym p (/C p ) of nonnegative divisors of degree p. 

Proof. Let M.00 be the set of divisors in Sym p (/C p ) for which (at least) one point 
is equal to Poo + - The image a P (Moo) oi Moo is given by 

a Po (Moo) = Ap^FooJ + « Po (Sym p - 1 (/C p )) C J(JC p ). (4.3) 

Since the (complex) dimension of Sym p_1 (K. p ) is p— 1, its image must be of measure 
zero by Sard's theorem (see, e.g., [U Sect. 3.6]). Similarly, let A^sp be the set of 
special divisors, then its image is given by 

ap (^ sp ) = ap (Sym p - 2 (A: p )), (4.4) 

assuming Pq to be a branch point. In particular, we conclude that cx Po (M sp ) C 
a Po (M.oo) and thus a P(j (A'ising) = fip (A / t O o) has measure zero, where 

X sing = Moo UM sp . (4.5) 

Hence, 

(J ("PoC^Ug) + «4p„_ ( p °o + )) (4.6) 
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is of measure zero as well. But this last set contains all initial divisors which will 
hit Poo + or become special at some n. We denote by Mo the inverse image of the 
complement of the set (14.61) under the Abel map, 

M = a p l UymP{K, p )\ |J (a Po (M sing ) + nA Paa _ (P^jA (4.7) 

Since Mo is of full measure, it is automatically dense in Sym p (/C p ). □ 

We briefly illustrate some aspects of this analysis in the special case p = 1 (i.e., 
the case case where (|3.1|) represents an elliptic Riemann surface) in more detail. 

Example 4.3. The case p = 1. 
In this case we have 

Fi(z,n) = z - /xi(n), 

G 2 (z,n) = i? 4 (Ai(«)) 1/2 + {z- b{n))F x {z,n), (4.8) 

3 

Ra{z) = Y[(z-E m ), 
and hence a quick calculations shows that 

\2 n _ /l„/„\2 



(4.9) 



G%(z, n) ~ Ri(z) — 4a(n) (z — fj,i(n))(z — H\{n + 1)) 

= (z - ^i(n))(4a(n) 2 z - 4a(nfb(n) + E), 
where 

E = \{Eo + E1-E2- E 3 )(E -E x +E 2 - E S )(E -E 1 -E 2 + E 3 ). (4.10) 
o 

Solving for fj,i(n+ 1), one obtains 

Ml (n + l) = 6(n)-^. (4.11) 

This shows that /ii(no + 1) — > 00, in fact, /ii(no + 1) = O(a(no) -2 ) as a(no) — > 
during an appropriate deformation of the parameters E m , m = 0,...,3. In 
particular, as a(no) — > 0, one thus infers b(no + 1) — > 00 during such a deformation 
since 



1 3 

6(n) = 5 £ S ™ - neZ, (4.12) 



ro=0 

specializing to p = 1 in the trace formula (|3.25|) . Next, we illustrate the set M^ 
in the case p = 1. (We recall that .M S p = and hence A4 s i ng = M.00 if p = 1-) By 
(|4.2[) one infers 

Ap^CM™)) = ^-P<x> + (Ai(^o)) + (n-n )A Poo+ (P oo _), n,n eZ. (4.13) 
We note that /fij G A^oo is equivalent to 

there is an n G Z such that /«i(n) = Poo + (or Poo_). (4-14) 
By (|4.13p . relation (|4.14|) is equivalent to 

J 4p oo+ (AiM)+A Poo+ (P oo _)Z = (modLO. (4.15) 
Thus, D/i^no) € Mo C /Ci if and only if 

J 4 Poo+ (Ai(no))+A Foo+ (P 00 _)Z^O (mod Li) (4.16) 
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or equivalently, if and only if 

A Poo _ (AiM) + A Poo _ (Poo + )% + (mod L{). 



(4.17) 



Next, we describe the stationary Toda algorithm. Since this is a somewhat 
lengthy affair, we will break it up into several steps. 
The Stationary (Complex) Toda Algorithm: 

We prescribe the following data: 
(i) The set 

{Emf^tl C C, E m ± E m > for to + to', to, m' = 0, . . . , 2p + 1 (4.18) 

for some fixed peN. Given {E m }^ , we introduce the function i?2 P +2 and the 
(nonsingular) hyperelliptic curve IC p as in (|3.1[) . 
(it) The nonspecial divisor 

%„ o) € SymP(/C p ), (4.19) 

where A( n o) is of the form 

Mno) = {Ai("o)> ■ ■ ■ , A P («o)} 

= {Al("o), ■ ■ ■ , Al(™o), ■ ■ ■ (Aa(no)) ■ ■ ■ )Ag(no)} 



(4.20) 



Pi ( n o) times 



P g (n )(«o) times 



with 



•>9(«o), 
(4.21) 

(4.22) 



A*("o) = (Mfc( n o)i2/(Afc("o))), Mfc(^o) 7^ Mfe'( n o) for ^ *') fe,fe' = l, 
and 

9("o) 

Pfc(n ) e N, = 1, . . . , g(rao), 2>fc(n ) = P- 

fc=i 

With {£ m }^ and 2? 

A(no) prescribed, we next introduce the following quantities 

(for zeC): 

q(n ) 

F p (z,n Q )= J] («-/i*(no)) Pfc(Bo) , 



fe=i 



? K)P,(no)-l (^(i2 2p+2 (C)V2)/^ 

r p -i(g,n ) = --F p (z,n ) > ^ 2^ — 



k=l 1=0 



i\(p k (n Q )-t-l)\ 



«("o) 



\ k' = l,k'^k / 



_^ 2p+l ?(«o) 

& (™o) = x J! ^ m ~~ X! Pk( n o)v>k(no), 



(4.23) 
(4.24) 

(4.25) 
(4.26) 



m=0 fc=l 

G p+ i(z,n ) = ~(z - b(n ))F p (z,n a ) +T p _i(z,n ). 
Here the sign of the square root in (|4.24p is chosen according to (|4.21[) . 

Afcfao) = (Mfc(^o),y(Afc(«o))) = (^k(n ), R 2p+2 {p.k(n )) 1/2 ) , k = 1, . . .,q{n Q ). 

(4.27) 

Next we record a series of facts: 
(I) By construction (cf. Lemma fB . 1 [) . 
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T ' _i(/Xfe(n ),no) 



= Gpli(fJ-k(n ),n ), 



(4.28) 



^ = 0, ... ,Pk(n ) - 1, fc = 1, ... ,g(no), 
(here the superscript (£) denotes £ derivatives w.r.t. z) and hence 

Afc("o) = (^k(n ), -G p+1 (nk(n ),n )), k = 1, . . . , g(n Q ). (4.29) 

(II) Since £>^( no ) is nonspecial by hypothesis, one concludes that 

Pfe(^o) > 2 implies R 2p+2 {^k(no)) / 0, fc = 1, . . . , q(n Q ). (4.30) 

(III) By (I) and (II) one computes 



dz e 



d e R 2p+2 (z) 



z=«(n ) a/ - 



z=M"o)' (4.31) 



z G C, I = 0, .. . ,Pk(n ) - 1, fe = 1, .. . ,g(n ). 

(IV) By (j4~26]) and ([4~3i"jl one infers that F p divides i? 2f >+2 - G 2 +1 . 

(V) By (j4~25]) and (|4~26]) one verifies that 

R 2p+2 (z) - G p+l (z,n a ) 2 = 0(z 2 ' } ). (4.32) 

z— >oc 

By (IV) and (|4.32|) we may write 

R2 P +2(z) - G p+1 (z,n ) 2 = F p (z,n Q )F p ^ r (z,n + 1), z G C, (4.33) 

for some r S {0, . . . , p}, where the polynomial F p - r has degree p — r. If in fact _Fb = 
0, then R 2p + 2 {z) = G p+ i(z,no) 2 would yield double zeros of R 2p + 2 , contradicting 
our basic hypothesis (|4.18ll . Thus we conclude that in the case r = p, Fq cannot 
vanish identically and hence we may break up (|4.33j) in the following manner 

i fT3 x V - G p+ i(z,n ) F p ^ r (z,n + 1) 

4>{P,n ) = — r = — — r, P = (z,y) e 1C P . 4.34) 

Fp(z,n Q ) y + G p+1 (z,n ) 

Next we decompose 

p—r 

Fp„ r (z,n + l) = C~[[(z-ii 3 (n Q + l)), z G C, (4.35) 
j'=i 

where (7 G C\ {0} and {^ij(no + l)}?=i C C (if r = p we follow the usual convention 

and replace the product in (|4.35|) by 1). By inspection of the local zeros and poles 
as well as the behavior near Poo ± of the function cf>( ■ , no), its divisor, (<^( • , no)) , is 
given by 

(<£(•, n )) = X ) P oo+ A(n +l) - p P»o_A(no)' ( 4 - 36 ) 

where 

A(n + 1) = {Al(«0 + 1),.- -jAp-r^O + l),Poo+, ■ • .,-Poo + }- (4.37) 

V v ' 

r times 

In particular, 

f/Vno+i) is a finite divisor if and only if r = 0. (4.38) 

We note that 

Oq Pa("0+1)) = "q (%«„)) -4p 00 _(^oo + ), (4.39) 



THE TODA HIERARCHY INITIAL VALUE PROBLEM 



19 



in accordance with (|4.2I) . 

(VI) Assuming that (14.32| is precisely of order z 2p as z — ► oo, that is, assuming 
r — in (|4.33|) . we rewrite (|4. 33[) in the more appropriate manner 

R 2p +2(z)~G p+1 (z,n ) 2 = -4a(n ) 2 F p {z, n )F p (z,n + 1), z G C, (4.40) 

where we introduced the coefficient a(no) 2 to make F p ( ■ , uq + 1) a monic polynomial 
of degree p. (We will later discuss conditions which indeed guarantee that r = 0, 
cf. (|4.38p and the discussion in step (XI) below.) By construction, F p ( ■ , no + 1) is 
then of the type 

g(no+l) s(no+l) 

F p (z,n + 1)= [] (z-Mnv + W k{no+1) , J2 P*("0 + 1)=P> f4 »x 
fc=i fc=i ^ ' ^ 

Ufc("o + 1) ^ "fc'("o + 1) for k^k' , k, k' = 1, . . . , g(n + 1), ^ e C, 
and we define 

Afe("o + 1) = (A»k(«o + l),G p +i(^fc(no + l),no)), fc = 1, . . . ,£?(no + 1). (4.42) 
Moreover, we introduce the divisor 

%„ Q+ D G SymP(/C p ) (4.43) 

by 

A(n + 1) = {/ti(n + 1), ■ • ., Ap(«o + 1)} 

= {Al( n + !),-■ - ,Al(«0 + 1),.. • ,A 5 (n +l)> ■ ■ -jAgCno+l)}- (4.44) 
v ^ ' . ' 

pi(«o+l) times P a (n + i)(no+l) times 

In particular, because of the definition (|4.42[) . f/i( no +i) is nonspecial and hence 

Pk{n + 1) > 2 implies ifep+a (Mk ("0 + 1)) ^ 0, fc = 1, . . .,g(n + 1). (4.45) 
Again we note that 

Oq (^a("o+i)) =%(%»)) -4p 00 _(-Poo + ), (4.46) 

in accordance with (|4.2p . 

(VII) Introducing 

_^ 2p+l «(rJo+l) 

6(n + 1) = - Em ~ p k {na + l)nk( n o + l), (4.47) 

m=0 k=l 

and interpolating G p +i( • , no) with i^( • , no + 1) rather than F p { ■ , no) yields 

G p+ i(z,n ) = -{z - b(n + l))F p (z,n + 1) -T p _i(z,n + 1), zeC, (4.48) 
where 

T p -i(z,n + 1) = F p (z,n + 1) 

fc fl(p fc (no + l)-i-l)l 1 J 

/ rfPkino+^-e-l ( 9("o + l) \ 

- ' pj-i (^-C)- 1 II (C-M fe '(no + l))- p *' ( " 0+1) 

\ fe'=l,fc'#fe // C=M"o + l) 
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Here the sign of the square root in (|4.49() is chosen in accordance with (|4.42|) , that 
is, 

/}fe(no + 1) = (nk(na + 1), y(Afe("o + 1))) 

= {fJtk(no + 1), Gp+i(fj,k{no + l),no)) 

= (« fe (n + l),i? 2p+2 (u fe (n + l)) 1/2 ), k = l,...,q{n + l). (4.50) 

(VIII) An explicit computation of a(n ) 2 then yields 

a(no) = 2 g feM^ 

9(«o) -, 

x 1] (i"feM-Mfe'(no))- pfc(no) + i (& (2) (n )-6(no) 2 ). (4.51) 
Here and in the following we abbreviate 

2p+l g(n) 

fe(2) w -E £ » - Ew( n )w(») J ( 4 - 52 ) 

m=0 fc=l 

for an appropriate range of n £ N. 

The result (|4.5ip is obtained as follows: One starts from the identity (14.40)) . 
inserts the expressions (I4.23|) and (|4.26[) for F p ( ■ , no) and G p+ \{ ■ , no), respectively, 
then inserts the explicit form (|4.24p of T p _i(- , no), and finally collects all terms 
of order z 2p as z — ► oo. An entirely elementary but fairly tedious calculation then 
produces (|4.51|) . 

In the special case q(no) — p, Pk(no) — 1, k= l,...,p, (|4.5ip and (|4.52|) reduce 
to (l3~2~oT) and (13T21) (for fc = 2). 

(IX) Introducing 

G p+1 (z,n Q + 1) = -(z - 6(n + l))F p («,no + 1) +T p _i(z,n + 1) (4.53) 
one then obtains 

G p+ i(z,na + 1) = -G p+ i(z,n ) - 2(z - 6(n + l))F p (z,n + 1). (4.54) 

(X) At this point one can iterate the procedure step by step to construct F p ( ■ ,ri), 
G P +i{ - , n), T p _i(< ,n), a(n), b(n), Ufc(n) , fc = 1, . . . , <?(n), etc., for n G [no, oo) DZ, 
subject to the following assumption (cf. ()4.38[) ) at each step: 

C/i(„+i) is a finite divisor (and hence a(n) ^ 0) for all n € [no, oo) n Z. (4.55) 

The formalism is symmetric with respect to n and can equally well be developed 
for n £ (— oo, no] H Z subject to the analogous assumption 

f A(n-i) i s a nm te divisor (and hence a(n) 7^ 0) for all n e (—00, no] (1 Z. (4.56) 

Indeed, one first interpolates Gp+i( ■ , no — 1) with the help of F p { ■ , no), then with 
F p (- ,n - 1), etc. 

Moreover, we once again remark for consistency reasons that 

Sqo^AO*)) =SLQ (PHn )) - (n-n )A Psxi _(P oo+ ), n e Z, (4.57) 
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in agreement with our starting point (14. 2|) . 

(XI) Choosing the initial data Vp,(n ) such that 

%„ ) 6 M), (4.58) 

where A4o C Sym p (/C p ) is the set of finite divisors introduced in Lemma T4.21 then 
guarantees that assumptions (|4. 55[) and (|4.56p are satisfied for all n G Z. 

(XII) Performing these iterations for all n G Z, one then arrives at the following 
set of equations for F p and G p +i after the following elementary manipulations: 
Utilizing 

G 2 p+1 - Aa 2 F p F+ = R 2p+2 = {G p+l f - A( a -) 2 F p F p , (4.59) 

and inserting 

G+ + i = -G p+1 - 2(z - b+)F+ (4.60) 
Gp+i - (G- +1 ) 2 - 4a 2 F p F+ + 4(a-) 2 F"F p = (4.61) 



into 



then yields 



2a 2 F+ - 2(a-) 2 F p - + (z - b)(G p+1 - G p+1 ) = 0. (4.62) 

Subtracting (|4.60|) from its shifted version G p +i = — G p+1 — 2(z — b)F p then also 
yields 

2(z ~ b+)F+ - 2(z - b)F p + G+ +1 - G p+1 = 0. (4.63) 

As discussed in Section^ (]4.62[) and (|4.63[) are equivalent to the stationary Lax and 
zero-curvature equations (|2.15[) and (|2.60[) . At this stage we have verified the basic 
hypotheses of Section [3] (i.e., Hypothesis 13.11 and the assumption that a,b satisfy 
the pth stationary Toda system (|2.20p ) and hence all results of Section [3] apply. 
Finally, we briefly summarize these considerations: 

Theorem 4.4. Suppose the set {-E m }m=o ^- ^ satisfies E m ^ E m > for m ^ ml , 
m, w! = 0, . . . , 2p+ 1, and introduce the function R 2p +2 and the hyperelliptic curve 
JC P as in f|3 . 1 [) - Choose a nonspecial divisor T>^ no ^ G M.o, where A4q G Sym p (/C p ) 
is the set of finite divisors introduced in Lemma \A.2\ Then the stationary (complex) 
Toda algorithm as outlined in steps (I)-(XII) produces solutions a,b of the pth 
stationary Toda system, 

s-TUa, b) = ( f P +1 ~ fp _ +1 ) = 0, p G N , (4.64) 
\9p+i ~9 p+ iJ 

satisfying (|3.2p and 



a(nf 



1 



«W (S(R 2p+2 (zfl 2 )jdz" 



\z=p, k (n) 



2 ^ (Pk(n) - 1)! 

?(«) , 
x II (Mfe(«)-^W)^ fc(n) + -(b^(n)-b(n) 2 ), (4.65) 
fe'=i, k'=ik 

^ 2p+l q(n) 

& (") = 2 ^ Em _ *52 Pk ( n ^ k ( n }> n eZ. (4.66) 

m=0 fe=l 

Moreover, Lemmas 13.2113.41 apply. 
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Remark 4.5. Suppose the hypotheses of the previous theorem are satisfied and 
a (no), b(n ), b(n + 1), F p (z,n ), F p (z,n + 1), G p+ i(z,n Q ), and G p+ i (z, n + 1) 
have already been computed using step (I)-(IX). Then, alternatively, one can use 

(a- fF- = a 2 F+ + 2~ 1 {z - b){G p+1 - G+ +1 ) + (z - b) 2 F p 

- (z - b+)(z - b)F+ , (4.67) 
G~ = 2((z - b+)F+ -{z- b)F p ) + G+ +1 (4.68) 

(derived from (|2.60|) ) to compute a(n), b(n), F p (z,n), G p +i(z,ri) for n < uq and 

a+F ++ = aF p - 2~ x (z - b)(G+ +1 - G p+1 ), (4.69) 
C++ =G- 2{{z - b++)F++ -{z- b+)F+) (4.70) 

to compute a(n — 1), b(n), F p (z,n), G p +\{z,n) for n > no + 1. 

Theta function representations of a and b can now be derived in complete analogy 
to the self-adjoint case. Since the final results are formally the same as in the self- 
adjoint case we just refer, for instance, to [6], [7], [9], [10], [TH Sect. 1.3], [18], [19], 
[20] . [2i] (cf. also the appendix written in [8]), [25], [30l Appendix, Sect. 9], (32] 
Sect. 9.2], [331 Sect. 4.5]. 

The stationary (complex) Toda algorithm as outlined in steps (I )-(XII ), start- 
ing from a nonspecial divisor T> ^( no ) G Mo, represents a solution of the inverse 
algebro-geometric spectral problem for generally non-self-adjoint Jacobi operators. 
While we do not assume periodicity (or even quasi-periodicity) , let alone real- 
valuedness of the coefficients of the underlying Jacobi operator, once can view this 
algorithm a continuation of the inverse periodic spectral problem started around 
1975 (in the self-adjoint context) by Kac and van Moerbeke [15], [16] and Flaschka 
[12], continued in the seminal papers by van Moerbeke [24], Date and Tanaka [7], 
and Dubrovin, Matveev, and Novikov [10] , and further developed by Krichever [15] . 
McKean [23], van Moerbeke and Mumford [25], Mumford [26], and others, in part 
in the more general quasi-periodic algebro-geometric case. 

We note that in general (i.e., unless one is, e.g., in the special periodic or self- 
adjoint case), Z?/i(n) will get arbitrarily close to Poo + since straight motions on the 
torus are generically dense (see e.g. [2] Sect. 51] or [17l Sects. 1.4, 1.5]). Thus, 
no uniform bound on the sequences a(n),b(n) exists as n varies in Z. In particu- 
lar, these complex- valued algebro-geometric solutions of some of the equations of 
the stationary Toda hierarchy, generally, will not be quasi-periodic (cf. the usual 
definition of quasi-periodic functions, e.g., in [311 p. 31]) with respect to n. For 
the special case of complex-valued and quasi-periodic Jacobi matrices where all 
quasi-periods are real- valued, we refer to [4] (cf. also [3]). 

5. Properties of Algebro-Geometric Solutions 
of the Time-Dependent Toda Hierarchy 

In this section we present a quick review of properties of algebro-geometric so- 
lutions of the time-dependent Toda hierarchy. Since this material is standard we 
omit all proofs and just refer to [5] (cf. also [H Sect. 1.4], [31 Chs. 12, 13]) for 
detailed presentations and an extensive list references to the literature. 

Throughout this section we will make the following assumption: 
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Hypothesis 5.1. Suppose that a,b satisfy 

a(-,t),b(-,t) £ C z , t £ R, a(n, •),&(»> •) G C X (R), neZ, 

a(n,t) ^ 0, (n,()eZxI ( ~ 5 ' 1 '' 

and assume that the hyperelliptic curve IC p , p G No, is nonsingular. 

In order to briefly analyze algebro-geometric solutions of the time-dependent 
Toda hierarchy we proceed as follows. Given p G No, consider a complex- valued 
solution a^°\b^ of the pth stationary Toda equation s-Tl p (et, 6) = 0, associated 
with IC P and a given set of summation constants {ci}e=i p C C. Next, let r G No; 
we intend to consider solutions a = a(t r ),b = b(t r ) of the rth Tl flow Tl r (a, b) — 
with a(t Q . r ) = a(°),6(t , r ) = 6 (0) for some io,r G R. To emphasize that the 
summation constants in the definitions of the stationary and the time-dependent 
Tl equations are independent of each other, we indicate this by adding a tilde on all 
the time-dependent quantities. Hence we shall employ the notation P2 r +2, 
F r , G r+ i, / s , g si c s , in order to distinguish them from P 2p +2, V p+ i, F p , G p+1 , 
fi, gi, ci, in the following. In addition, we will follow a more elaborate notation 
inspired by Hirota's r-function approach and indicate the individual rth Tl flow 
by a separate time variable t r G HL More precisely, we will review properties of 
solutions a, b of the time-dependent algebro-geometric initial value problem 

Th(a, 6 ) = h -<^ a 'S -f^t^) = 0, 

V K+9p+i\aM -g p+1 {a,b) J ( 5 2 ) 

M)U ,^ (0 ^ (0) ), 

S " Tlp(a ' 6 } -{ 9p+1 (a(°),6(°))- 9p+1 (a(°),#)) )-° (5 ' 3) 

for some to.r £K,p,r£ No, where a = a(n, t r ), b — b(n, t r ) satisfy (|5.ip and a fixed 
curve K,p is associated with the stationary solutions a^°\b^ in (|5.3p . In terms of 
Lax pairs this amounts to solving 

~ [?2r+2(t r ),L(t r )} =0, t r £ M, (5.4) 

[P2H-2(t0,r),£(t0,r)]= - ( 5 ' 5 ) 

Anticipating that the Tl flows are isospectral deformations of L(to >r ), we are going 
a step further replacing (|5.5p by 

[P2 P+ 2(tr),L{t r )]=0, t r eR. (5.6) 

This then implies 

2p+l 

P2p+2(trf = R2p+2(L{t r )) = Y[ (L(t r ) - E m ), t r G R. (5.7) 

Actually, instead of working with (|5.4p , (|5.5p , and (|5.6p , one can equivalently take 
the zero-curvature equations (|2.63p as one's point of departure, that is, one can also 
start from 

U tr + UV r+1 - V+ +1 U = 0, (5.8) 
UV p+1 -V+ +1 U = 0, (5.9) 
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where (cf. f2T23]t . (l2~24l) . (l2~58)l . (l2~59l) s ) 

1 



t/(z) 



—a /a (z — b)/a 



Vh-i(«) = 



-2a-F r {z) 2{z - b)F r {z) + G r+1 {z) 



and 



^) = Z)/p-^ = II( js -^)' /o = 1 ' ( 5 - n ) 

p 

Gp + i(z) = -z p+1 + ^.g p _^ + / p+1 , flo = -ci > (5.12) 



F r (z) = J2fr-sZ s , /o = l, (5.13) 

r 

G r +i(z) = -z r+1 + ^.g r _ s z s + / r+ i, ff = -ci, (5.14) 

s=0 

for fixed p,r 6 No- Here /g, / s , g^, and </ a , ^ = 0, . . . ,p, s = 0, . . . , r, are defined 
as in (|2.4[) - (|2.6[) with appropriate sets of summation constants eg, £ £ N, and Ck, 
k e N. Explicitly, (EU) and (EU) are equivalent to (cf. (j2"35)l . J33BJ), (|232]1 . (^55]) ) 

a tr - -a(2(* - b+)F+ + G+ +1 + G r+1 ), (5.15) 

b tr = 2((z - b) 2 F r + {z- b)G r+l + a 2 F+ - {a-) 2 F~), (5.16) 

= 2(z-6+)F+ + G+ +1 + G p+1 , (5.17) 

= (z - 6) 2 F p + (z - 6)G p+ i + a 2 F+ - (a") 2 ^, (5.18) 

respectively. In particular, (|2.34|) holds in the present independent setting, that is, 

G 2 +1 - Aa 2 F p F+ = R 2p+2 . (5.19) 

As in (13.61) one introduces 



fLj(n,t r ) = {/j,j(n,t r ), -G p+ iQij(n,t r ),n,t r )) £ fC p , j = 1, . . . ,p, (n,i r ) e Z x K, 

(5.20) 

p,f(n,t r ) = (n,t r ),G p+1 (fxj (n,t r ),n,t r )) £ fC p , j = 1,... ,p, (n,t r ) e Z x R, 

(5.21) 

and notes that the regularity assumptions (|5.ip onai imply continuity of /ij with 
respect to t r 6 R (away from collisions of zeros, jiij are of course C°°). 

In analogy to (|3.7[) , (|3.8p . one defines the meromorphic function tfi(-,n,t r ) on 

2a(n, t r )b p {z, n, t r ) 
_ -2a(n, U)F p (z, n+l, t r ) ^ 



y + G p+ i{z,n, t r ) 
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P{z,y) e Kp, (n,t r ) e Z x R, 
with divisor (cj>( ■ , n, t r )) of 4>( ■ , n, t r ) given by 

(4>(- ,n,t r )) =% + g„+i, tr )-I'p co j(„,t r ), (5-24) 

using (pTTTj) and (|QD| . 

The time-dependent Baker- Akhiezer function ij)(P, n, no, t r , to, r ), meromorphic 
on IC p \ {Pqc + 7 Poo- }i is then defined in terms of <j) by 

^p{P,n, n ,t r ,to, r ) 

= exp ^ ds(2a(n ,s)F r (z,rio,s)(f>(P,no,s) + G r+ i(z,n ,s))^j 

{rim^no ^( P ' m ' *»■) for n > n + 1, 
1 for n = n , (5.25) 

IC^ 1 <MP m, U)- 1 for n < no - 1, 
P&)C P \{Poo ± }, (n,no,t r ,MeZ 2 xM 2 . 
For subsequent purposes we also introduce the following Baker-Akhiezer vector, 

*(P,n,fH, > t r ,V)=^/7p'^^^ )> ). 
v V i>{P,n,n ,t r ,ta,r) J (5.26) 

i'GX^\{P 00± }, (n,n ,t r ,t ,r) e Z 2 xM 2 . 

Basic properties of 0, ?/>, and \& are summarized in the following lemma. 

Lemma 5.2. Assume Hypothesis 15.11 and suppose that a,b satisfy (I5.15P (|5.18|) . 

In addition, let P — (z,y) € /C p \ {Poo ± } 7 (n,n ,t r ,to yr ) G Z 2 x K 2 7 and r 6 No- 
TTien satisfies 

a<t>{P) + a~ (4>- (P))- 1 = z - b, (5.27) 
J> tr (P) = -2a(P r (z)0(F) 2 + F+(z)) + 2(z - b + )F+ (z)<j>(P) 

+ (G+ +1 (z)-G r+1 (z))4>(P), (5.28) 

m4>(P*) = fr^, (5-29) 
#P)-0(P*) = Jjj£L, (5.30) 

Moreover, ip and $ satisfy 

(L - z(P))ip(P) — 0, (P 2p+2 - y(P)MP) = 0, (5.32) 

^ r (P)=P 2r + 2 V(P) ( 5 - 33 ) 
= 2aP r (z)^+(P) + G r+1 (z)^(P), (5.34) 

*+(P) = C/(z)*(P), y*(P) = F p+ i*(P), (5.35) 

* tr (P)=V P+ i(«)*(P), (5.36) 

■0(P, n, n , i r , t 0r )ip(P*,n, n Q ,t r , t 0r ) = Fp( z > n ' t r) (5.37) 

F p (z,n ,t ,r) 
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a(n, t r ) (ip(P, n, n , t r , to, r )ip(P* , n + 1, n , U, to, r ) 

+ ip(P* ,n,n ,t r ,t 0<r )ip(P,n + l,n ,t r ,t Q . r )) = _ Gp+i(z,n,t r ) (5.38) 

F p (z,n ,t 0tr ) 

W(j>(P, ■ , n , t r , t , r ), 1>{P*, ■ , no, U, to,r)) = ~ F / (P \ c - (5-39) 

In complete analogy to the case of stationary trace formulas one obtains trace 
formulas in the time-dependent setting (cf. the abbreviation (|3.24[l ). 

Lemma 5.3. Assume Hvvothesis 15.11 and suppose that a,b satisfy (I5.15[) f)5.16[) . 

Then, 

2p+l p 

6= 2^ £m_ ^- (5 - 40) 

m— j— 1 

In addition, if for all (n, t r ) £ Z x K, (J,j(n, t r ) ^ ^k(n, t r ) for j ^= k, j,k = 1, . . . ,p, 
then, 

« 2 = ^X>0%) 11^ - + J(& (2) - 6 2 )- (5.41) 

j=l k=l 

For completeness we next mention the Dubrovin equations for the time variation 
of the Dirichlet eigenvalues of the Toda lattice. 

Lemma 5.4. Assume Hvvothesis 15.11 and suppose that a,b satisfy (|5.15|) - (l5.16p . 

In addition, assume that the zeros /ij(n,t r ), j = l,...,p, of F p ( - ,n,t r ) remain 
distinct for all (n,t r ) 6ZxR. Then, 

d . . ~ , , . . y(fij(n,tr)) 

—Hj{n,t r ) = -2F r (fj,j{n,t r ),n,t r ) p 

dt r ll|=i(MMr) -mn,U)) (542) 

j = 1, . . . ,p, (n,t r ) e Z x K. 

When attempting to solve the Dubrovin system (|5.42[) . it must be augmented 
with appropriate divisors T>^ no to \ 6 Sym p /C p as initial conditions. 

For the independence of F p and G p +i one obtains the following result. 

Lemma 5.5. Assume Hypothesis 15.11 and suppose that a,b satisfy (|5.15j) - (|5.16p . 

In addition, let (z,n,t r ) ECxZxR. Then, 

F p ,t r = 2(F p G r +i — Gp+iF r ), (5.43) 

G p+lttr = Aa 2 (F p F+ - F+F r ). (5.44) 
In particular, (|5.43l) and (|5.44ll are equivalent to 

V P+ i,t r = [V r+ i,Vp +1 ]. (5.45) 

It will be shown in Section[6]that Lcmma [5.5l in conjunction with the fundamental 
identity (|5 . 19|) yields a first-order system of differential equations for fi, gi, £ = 
l,...,p, that serves as a pertinent substitute for the Dubrovin equations (|5.42|) 
even (in fact, especially) when some of the /Xj coincide. 

As in the case of trace formulas, also Lemma l3.4l on nonspecial Dirichlet divisors 
T>jx and the linearization property of the Abel map when applied to extends 
to the present time-dependent setting. For the latter fact we need to introduce a 
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particular differential of the second kind, fl r , defined as follows. Let ui p ± q be the 
normalized Abelian differential of the second kind (i.e., with vanishing a-periods) 
with a single pole at Poc± of the form 

w p1 ± ,«= (C" 2 ^ + °( 1 )) d CnearP co± , <; £ Mo- (5.46) 
Given the summation constants c\, . . . , c r in F r (cf. (|5.13[) ). we then define 



S]f)^( ? + l)c r _,(^ +i9 -^ J, 20 = 1. (5.47) 

9=0 



(2) 

Since the differentials uj,^ q were supposed to be normalized we have 

[ f2^=0, j = l,...,p. (5. 

Jcij 

Moreover, writing 

^ = ( £ ^,m(Poo ± )C ro ) d( = ±( jr d,, m (Poo + )C m ) d( near P oc± 

^ m=0 ' ^ m=0 



(5.49) 



relation (|A.20|) yields for the vector of ^-periods U^ r ^ of Sl^ , 

i(2) 



U_y = (U%,...,U<$), (5-50) 



™ 2 



1 /■ ~ r 

- / ft( 2 > =2Vc r _,d j ,,(P 00+ ) ) j = l,...,p. (5.51) 

The time-dependent analog of Lemma EOl then reads as follows. 

Lemma 5.6. Assume Hypothesis 15.11 and suppose that a,b satisfy ([5~T5l) (f5~16|) . 

LefD^, fj, — {/xi, . . . , G Sym p (/C p ), &e £/ie Dirichlet divisor of degree p associated 
with a,b defined according to (|5.20|1 . that is, 

(ij{n,t r ) = (iu,j(n,t r ), -G p+ i(^(n,i r ),n, t r )) e /C p , j = 1, . . . ,p, (n,i r ) e Z x R. 

(5.52) 

T/ien 'Dp,(n,t T ) * s nonspecial for all (n, t r ) € Zxl. Moreover, the Abel map linearizes 
the auxiliary divisor T)^ in the sense that 

~(2) 

"Qo (%",*>•)) = SQo^ACno.to.r)) - (« - n o)i.P 00 _ (-Poo+J ~ (*r ~ t ,r)U r , (5.53) 

~(2) 

where Qq S /C p is a given base point and U_ r is the vector of b-periods of the 
differential of the second kind Qf- introduced in (|5.51|) . 

In addition, if a,b e L°°(Z x R), i/ien i/iere exists a constant C M > smc/i i/iai 

|A*i(n,«r)| < C MJ j = l,..-,p, (n,t P ) GZxR. (5.54) 
Proof. We will prove that 

V(P, n, n , t r , t , r ) =C(n, t r ) f ' ^ r)) ., (5.55) 
x exp f (n- n ) / ^ P + (t r - t , r ) / ] , 
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where 

z(P, n,t r ) = A Qo (P) - a Qo (%„ ;tr) ) + S Qo . (5.56) 

By Lemma 13.4 of |3j2 it suffices to show that the essential singularities at Poo ± 
are equal. That is, by (|5.25j) we need to show that 

ip(P, n ,no,t r ,t , r ) 

= exp ^jf ds(2a(n ,s)F r (z,n ,t r )(f>(P,n ,to :r ) + G r +\(z, n , s))^ 

= exp ( ±(t r -to, r )'f2cr-kC k ~ 1 +0(1) J for P near P^. (5.57) 
V fe=0 / 

Using ([5~22]) and ([Q3]) one obtains 

,/ D . , n / F p (z,n Q ,t r ) \ 1/2 / /" tr F r (z,n ,s)\ 

i/j{P,no,no,t r ,t ,r) = TT7 7 — r ) exp W ds — (5.58) 

\F p (z,n ,to,r)J \ Jt 0r Fp(z,n ,s) 



and the desired asymptotics follow from Theorem IC.1[ which tells us that 

F k (z,n , s) = ±C k ~ 1 + 0(1) for P near P OQ± , (5.59) 



F p (z, n , s) 

together with (|2^5|) . □ 
Again the analog of Remark 13.51 applies in the present time-dependent context. 

6. The Algebro- Geometric Toda Hierarchy Initial Value Problem 

In this section we consider the algebro-geometric Toda hierarchy initial value 
problem (|5.2[) . (|5.3[) with complex- valued initial data. For a generic set of initial 
data we will prove unique solvability of the initial value problem globally in time. 

While it is naturally in the special self-adjoint case to base the solution of the 
algebro-geometric initial value problem on the Dubrovin equations (|5.42p (and the 
trace formula (|5.40|) for b and formula (|5.41[) for a 2 ), this strategy meets with 
difficulties in the non-self-adjoint case as Dirichlet eigenvalues p,j may now collide 
on ICp and hence the denominator of (|5.42|) can blow up. Hence, we will develop 
an alternative strategy based on the use of elementary symmetric functions of the 
variables {(J'j}j=i,...,p in this section, which can accommodate collisions of flj. In 
short, our strategy will be as follows: 

(i) Replace the first-order autonomous Dubrovin system (|5.42[) of differential 
equations in t r for the Dirichlet eigenvalues fij(n,t r ), j — 1, . . . ,p, augmented by 
appropriate initial conditions, by the first-order autonomous system (|6.27l) . (|6.28[) 
for the coefficients fj, j — 1, . . . ,p, gj, j = 1, . . . ,p— 1, and g p + f p+ i with respect 
to t r . (We note that fj, j = l,...,p, are symmetric functions of . . . , ji p .) 
Solve this first-order autonomous system in some time interval (to,r — To, to rT + Tq) 
under appropriate initial conditions at (no, fo, r ) derived from an initial (nonspecial) 
Dirichlet divisor £>£(„ Mo r ). 

(ii) Use the stationary algorithm derived in Section 2] to extend the solution of 
step (i) from {no} X (to,r — Tq, to, r + Tq) to Z x (i , r — Tq, to, r + Tq) (cf. Lemma 

E2J. 

(Hi) Prove consistency of this approach, that is, show that the discrete algo- 
rithm of Section [4] is compatible with the time-dependent Lax and zero-curvature 
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equations in the sense that first solving the autonomous system ()6.27|) . (|6.28|) and 
then applying the discrete algorithm, or first applying the discrete algorithm and 
then solving the autonomous system (|6.27|) , (|6.28p yields the same result whenever 
the same endpoint (n, t r ) is reached (cf. the discussion following Lemma 16.21 and 
the subsequent Lemma 16.31 and Theorem I6.4[) . 

(iv) Prove that there is a dense set of initial conditions of full measure for which 
this strategy yields global solutions of the algebro-geometric Toda hierarchy initial 
value problem (cf. Lemma 16.51 and Theorem 16. 6[) . 

To set up this formalism we need some preparations. From the outset we make 
the following assumption. 

Hypothesis 6.1. Suppose that 

aieC z and a(n) / for all n e Z, (6.1) 

and assume that a,b satisfy the pth stationary Toda system (|2.20p . In addition, 
suppose that the hyperelliptic curve JC P in (|3.ip is nonsingular. 

Assuming Hypothesis 16. 1[ we consider the polynomials F p , G p +\, F r , and G,+i 
given by (|5.11|) -(|5.14|) for nx ed p, r G N . Here ft, f s , g e , and g s , t = 0, . . . ,p, s = 
0, . . . , r, are defined as in (|2.4[) - (|2.6[) with appropriate sets of summation constants. 

Our aim will be to find an autonomous first-order system of ordinary differential 
equations with respect to t r of ft and ge rather than fij. Indeed, we will take the 
coupled system of differential equations (|5.43|) . (|5.44p . properly rewritten next, as 
our point of departure. In order to turn (|5.43p . (|5.44p into a system of first-order 
ordinary differential equations for fi and gi, we first need to eliminate f s , g s , 
and /J 1 " in terms of fe and gi as follows. 

Using (|2ll), (|2~23| . (|2~25l) . and Theorem [CT] one infers 

r r 

^W = ^/r-^ S =^£r-s^(z), (6.2) 
s=0 s=0 

I lAp 

F e (z) = J2fe-kz\ /o = l, h = Y,tt-k(E)f k , leMo, (6.3) 

A:=0 k=0 

where m An = min{m,n} and ci{E) has been introduced in (|C4[) . Hence one 
obtains 

/o = l, / s = ^i, s (/i,... J /p), s=l,...,r, (6.4) 

where J-i, s , s = 1, . . . , r, are polynomials in p variables. 

Next, using (pEHj) . ([2~2%]) . and Theorem [CT] one concludes 

r r+1 

G r+1 (z) = -z r+1 + ^.g r _ s z s + / r+ i = y^ y c r+1 - s G s (z), 

s=0 s=l 

G (z) = G (z)| ci=0 = 0, G^z) = G x {z) = -z-b, 

i 

Ge+i{z) = G e+1 {z)\ ck=ok ^^ l = -z l+1 +Y,k-kz k ^N, (6.5) 

fc=0 

lAp 

9o = 0, ge = ^ct- k (E)(gk + f P+ iS p , k )-ce +1 {E), leN. (6.6) 

fc=0 



30 



F. GESZTESY, H. HOLDEN, AND G. TESCHL 



Hence one concludes 

9o = -ci, g s = J r 2,s(fi,---,f P ,9i,---,9p-u(gp + fp+i)), s = l,...,r, (6.7) 

where ^si s = 1, . . . ,r, are polynomials in 2p variables. We also recall (cf. (|2.18|) ) 
that /p+i is a lattice constant, that is, 

/p+i = / P +i- (6-8) 
Next we invoke the fundamental identity (|2.34|) in the form 

-4a 2 F/ = J?2p+2 ~^ +1 . (6.9) 

While (|6.9|) at this point only holds in the stationary context, we will use it later 
on also in the independent context and verify after the time-dependent solutions 
of (|5.2|) , (|5.3[) have been obtained that (|6.9J) indeed is valid for all (n,t r ) elxl. 
A comparison of powers of z in (|6.9[) then yields 

4a // = T 3 j(fi, ...,f p ,gi,.. .,g p -i, (g p + f P +i)), £=l,...,p, 

where J-^.e, i = 1, . . . ,p, are polynomials in 2p variables. 

Finally, combining ([5^]) . fO|) . fO|) . and (jCTT))) . one obtains 

(611) 

4a 2 /+ = ^4, s (/i, ■ ■ .,fp,gi,...,g p -i, (g p + / P +i)), s = 1, . . . ,r, 

where J-4. s , s = 1, . . . , 3, are polynomials in 2p variables. 

We emphasize that also the Dubrovin equations (|5.42p require an analogous 
rewriting of F r in terms of (symmetric functions of) fj,j in order to represent a 
first-order system of differential equations for fj,j, j = 1 , . . . , p. 

Next, we make the transition to the algebro-geometric initial value problem 
(|5.2|) . (15. 3p . We introduce a deformation (time) parameter t r G R in a = a(t r ) 
and b = b(t r ), and hence obtain independent quantities fi — fe(t r ), gg — gt(t r ), 
F p (z) = F p (z,t r ), G p +\{z) = G p +i(z,i r ), etc. At a fixed initial time io,r G R we 
require that 

(o,6)|t r=tG , r - (a (0) ,& (0) ), (6.12) 

where 

■ i to tr ), b^ — b( - ,to,r) satisfy the pth stationary Toda equation 
(|2.20p as in (|6.ip - (|6.1ip . As discussed in Section [4] in order to guarantee that 
the stationary solutions (|6.12p can be constructed for all n € 1 one starts from a 
particular divisor 

%«„,*„,,.) e Mo (6.13) 

where jj,(no,to, r ) is of the form 

A("0,*0,r) = {Al(™0,*0,r, ■ ■ ■ , Ap(™0, *0,r)} (6.14) 
= {Al(^0,i0,r), • ■ ■ , /*l(rao,*0,r), • ■ ■ , Ag(n ,t r )(^0; ^O.r), ■ ■ • , Ag(n ,to r) ( n 0; *0,r)} 

v v ' : v . ' 

Pi(n ,to,f) times P 9 (n. ,t 0ir )(no,io,r) times 

with 

Afc("0,*0,r) = (^k(no,to,r),y(^k(no,t ,r))), 

Mfc( n o,*o,r) / Mfc'( n o,*o,r) for fc / fc', fc,fc' = l,...,g(no,to,r) 5 
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and 

q(no,to,r) 

Pk{noM,r) G N, k = l,...,q(no,to,r), Pk(no,to,r) = P- (6.16) 

fc=i 

Next we recall 

p g(no,to,tO 
F p (z,n ,t ,r) = ^2fp-e{no,t ,r)z e = J] (2 - fik(n , *o,r)) Pfe(no '*°''' ) , 



£=0 fc=l 



(6.17) 



T p -i(z,Tlo,to,r) — —Fp(z,nQ,to, r ) 



X 



ti h> t\(p k {n Q ,t Q , r )-t-l)\ 



I 4Pk{no,to,r)-i-l I 

x . - . , - (z - C)" 1 (6.18) 

9(no,to,r) \ 

k'=l,k'jtk J J C=M*("o,to,r) 

2p+l «(no,to,-r) 

H«o I to,r) = 2l] £m_ X! Pfc( n 0,*0,r)/-tfc(™0,*0,r), (6.19) 

m=0 fe=l 



G p+ i(z,n ,to,r) = -z p+1 + } j g p -i(no,to t r)z e + f p+ i(t , r ), 



e=o 



= — (z — b(no, t ,r))F p (z, n n , to, r ) + T p -i(z, no, t , r )- (6.20) 
Here the sign of the square root in (|6.18[) is chosen as usual by 

Afc( n 0»*0,r) = (w(BO,to,r).#t(lO.V))) 

= (^fc(«o,io,r),i?2 P +2(A*fe(^o,io,r)) 1/2 ) (6-21) 
= (i«fc(no,to,r)) -Gp + i(/Ufc(no,io,r)>rco,io,r))> fc = 1, . . . , q(nn,t ^ r ). 

By (|6.17l) one concludes that (|6.14[) uniquely determines F p (z, no, to, r ) and hence 

fi(no, to <r ), f p (n , *o, r )- (6-22) 

By f|6.18[) — (|6.22[) one concludes that also Gp+i(z, no, to,r) and hence 

ffl(no,*0,r)) • • ■>5p-l("0»*0,r),5p('^0)*0,r) + /p+l(*0,r) (6.23) 

are uniquely determined by the initial divisor T>^^ no tg \ in (|6.13|) . 

Summing up the discussion in (|6.2[) ~ (|6.23p . we can transform the differential 
equations 

Fp,t r { z > n o,U) = 2(F p (z,n ,t r )G r+1 (z,n ,t r ) 

- G p+ i{z,n Q ,t r )F r (z,n Q ,t r )), (6.24) 
Gp+i,t r [z, n , t r ) = 4a(n , t r ) 2 (F p (z, n , t r )F+(z, n , t r ) 

- F+(z,n ,t r )F r (z,n ,t r )) (6.25) 
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subject to the constraint 

a 2^+/ . \ R2 P +2(z) - G p+ i(z, n ,t r )' 2 , r .>,-, 

-4a FJ(z,no,t r ) = — * — — ^— , 6.26) 

F p {z,n ,t r ) 



and associated with an initial divisor X>/}(„ 0) t 0r ) in (|6.13p into the following au- 
tonomous first-order system of ordinary differential equations (for fixed n — n ), 

fj,t r = Fj(fi, • • • j fp, 9i, ■ ■ ■ : 9p-i,9p + fp+i), 3 = 1) • • • ,P, 
9j,t r = dj(fx, ...,f p ,gi,.. .,g p -i,g p + f p +i), j = 1, . • • ,p - 1, (6.27) 
{9p + fp+i)t r = G P (fi, ■ ■ ■ j Jp,9ii ■ ■ • 5 9 P -i,9p + /p+i) 
with initial condition 

fj(no,t 0tr ), j = l,...,p, 

9j(no,to, r ), i = l,---,p-l, (6-28) 

g P (no,to,r) + /p+l(to,r)) 

where J-'j, Qj, j = 1, . . . ,p, are polynomials in 2p variables. As just discussed, the 
initial conditions (|6.28|) are uniquely determined by the initial divisor X>n(„ 0i t r ) in 
(16331) . 

Being autonomous with polynomial right-hand sides, there exists a To > 0, such 
that the first-order initial value problem (|6.27p . (|6.28[) has a unique solution 

fj = fj( n 0,tr), j = l,...,p, 
9j = 9j(no,t r ), j = l,...,p-l, 
9 P + fp+i = 9p{no, U) + f p +i(t r ) 
for all t r G (to, r - T , to, r + To) 

(cf., e.g., [351 Sect. III. 10]). Given the solution (|6.29[) , we next introduce the fol- 
lowing quantities (where t r G (to, r — 2o,to, r + To)). 

p i(n ,t r ) 
F p (z,n ,t r )=Y,f P -e(ri Q ,t r )z e = J] (z - ,i fe (n , f r )f*("<»H (6.30) 

1=0 k=l 

T p -i(z,n ,t r ) = -F p (z,n ,t r ) 



x 



C = Mfc(«0,*r) 



fc=l 



^!(p fc (n ,t r )-f-l)! 



/ (Ph(n ,U)-i-l ( 



q(n ,t r ) > 

x n (c-Mfe'K,^))^'^ ^ 

k' = l,k'^k / 
2p+l q(n ,t r ) 



C=Mfc("0,*r) 



b(n ,t r ) = ^y^E m - p k (n ,t r )fi k (n ,t r ), (6.32) 



2 

m=0 fe=l 
P 



G p+ i(z,n ,t r ) = + ^ g p ^i{n , t r )z e + f p+ i(t r 



e=o 



-{z - b(n a ,t r ))F p (z,n ,t r ) + T p ^i(z, n Q , t r ). (6.33) 
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In particular, this leads to the divisor 

%„„,M £ Sym^/Cp) (6.34) 
and the sign of the square root in (16 . 3 1 1) is chosen as usual by 

ftk(n ,t r ) = (nk(no,t r ),-Gp +1 ((j, k (no,t r ),no,t r )) 

, (6.35) 
= (nk{no 7 t r ) 7 R 2 p + 2(Hk(no,t r )) 1 ), k = 1, . . . ,q(n ,t r ) 

and 

jj,{n ,t r ) = {fj,i(n ,t r , . . . ,{J, p {n Q ,t r )} 

= { /ii(n ,t r ), . . .,m(n ,t r ) , . . . ,n q ( no j r )(n ,t r ), . . . , H q ( no ,t r )(no,t r )} (6.36) 

p\(n .t T ) times P 9 (n ,t T ) (™0,tr) times 

with 

u fc (n ,i r ) / Hk'(no,t r ) for fc ^ fc', fc, fc' = 1, ... ,q(n ,t r ) (6.37) 

and 

g(n ,t r ) 

Pk{no,U) e N, fc = 1, . . . ,q(n ,U), pk(n ,t r ) = p. (6.38) 

By construction (cf. (|6 .35|) ) . the divisor 'D^i^^ is nonspecial for all t r 6 (io,r — 
2o, io,r + 7b). 

In exactly the same manner as in (|4.28p ~ l|4.3ip one then infers that F p ( ■ ,n ,t r ) 
divides R 2p +2 ~ G 2 +1 (since t r is just a fixed additional parameter). Moreover, 
arguing as in (|4.32|) - (|4.38j) we now assume that the polynomial 

R2 P+ 2{z)-G p+1 {z,n ,t r ) 2 = 0(z 2 P) (6.39) 



z — >oo 



is precisely of maximal order 2p for all t r G (io.r — 7b, £o,r + 7b). One then obtains 
R2 P +2(z) - G p+ i(z,n Q ,t r ) 2 = -4a(n ,t r ) 2 i ;l p(z,no,tr)^p(-Z,«o + Mr), 

(6.40) 

(z, i r )eCx (i 0: r - Tb,t ,r + T ), 

where we introduced the coefficient a(no,t r ) 2 to make F p ( - , no + l,t r ) a monic 
polynomial of degree p. As in Section 2J the assumption that the polynomial 
Fp{' ,no + l,t r ) is precisely of order p is implied by the hypothesis that 

%no,M e Mo for all i r G (to,r - T , t , r + T ), (6.41) 

a point we will revisit later (cf. Lemma 16 . 5|) . Given (|6.40p . we obtain consistency 
with (|6.9| for n = hq and t r G (io.r — 7b, £o,r + 7b). 

The explicit formula for a(no, t r ) 2 then reads (for t r € (to,r — 7b,io,r + 7b)) 



a(n ,i r ) =- > -— — — t^t 

2 {pk(n ,t r ) - ly. 



z=fj, k (n ,t r ) 



fe=l 

g(n, 0l *T-) 

x J] (n k (n ,t r )- n k ,(n ,t r ))-P»( n °^ (6.42) 



+ j(b^(n a ,t r )-b(n ,t r ) 2 ). 
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Here and in the following we use the abbreviation 

^ 2p+l q(n,t r ) 

b (2 \n,t r ) = g XX - Pk(n,t r )fi k (n,t r ) 2 (6.43) 

m=0 k=l 

for appropriate ranges of (n,t r ) £ N x M. 

With (I6.30p - (I6.43I) in place, we can now apply the stationary formalism as sum- 
marized in Theorem 14.41 subject to the additional hypothesis (|6.41j) . for each fixed 
t r G (to,r ~ To,to. r + To). This yields, in particular, the quantities 

F p , G p+ i, a, b, and ft for (n,t r ) € Z x (t , r - T ,t , r + To), (6.44) 

which arc of the form (|6.30|) — (|6.43|) . replacing the fixed no 6 Z by an arbitrary 
rt 6 Z. In addition, one has the following fundamental identities (cf. (I4.54| . (I4.59|) . 
()4.62|) . and (|4.63[) ) . which we summarize in the following result. 

Lemma 6.2. Assume Hypothesis 16.11 and condition (|6.41|) . Then the following 
relations are valid: 

Rip+i — Gp +1 + 4a 2 FpFp = 0, 
2{z - b+)F+ + G+ +1 + G p+1 = 0, 
2a 2 F+ - 2(a-) 2 F- + (z - b)(G p+1 - G~ +1 ) = 0, 
2{z - b+)F+ - 2{z - b)F p + G+ +1 - G~ +l = 
on C x Z x (t 0>r -T ,t 0>r + T ) 

and hence the stationary part, (|5.9[) . o/ the algebro- geometric initial value problem 
holds 

UV P+1 - V p + +1 U = on C x Z x (< , r - T ,i 0)r + T ). (6.49) 
In particular, Lemmas I3.2H3.4I apply. 

Lemma 16.21 now raises the following important consistency issue: On one hand, 
one can solve the initial value problem (|6.27p . (|6.28|) at n — no in some interval 
tr G (to,r— To, to >r +To), and then extend the quantities F p , G v +\ to all C xZx (to,r — 
To,to, r + To) using the stationary algorithm summarized in Theorem 14.41 as just 
recorded in Lemma 16.21 On the other hand, one can solve the initial value problem 
(|6.27p . (|6.28p at n = m, n\ ^ n , in some interval t r £ (t . r — Ti,to, r + T\) with 
the initial condition obtained by applying the discrete algorithm to the quantities 
F p ,G p +i starting at (no,to, r ) and ending at (ni,to, r ). Consistency then requires 
that the two approaches yield the same result at n = n\ for t r in some open 
neighborhood of to, r - 

Equivalcntly, and pictorially speaking, envisage a vertical i r -axis and a horizontal 
n-axis. Then, consistency demands that first solving the initial value problem 
(|6.27p . (|6.28p at n = no in some i,.-interval around to. r and using the stationary 
algorithm to extend F p ,G p +\ horizontally to n = n\ and the same t r -interval 
around to. r , or first applying the stationary algorithm starting at (no, £o,r) to extend 
F p , Gp+i horizontally to (ni,to, r ) and then solving the initial value problem (|6.27p . 
(|6.28p at n = n\ in some i r -interval around to, r should produce the same result at 
n — n\ in a sufficiently small open t r interval around to,r- 

To settle this consistency issue, we will prove the following result. To this end 
we find it convenient to replace the initial value problem (|6.27p . (|6.28p by the 



(6.45) 
(6.46) 
(6.47) 
(6.48) 
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original ^-dependent zero-curvature equation (|5.8|) . Ut r + UV r +i — V^-JJ = on 
CxZx(f , r -To,V+T ). 

Lemma 6.3. Assume Hvvothesis Vo.W and condition (|6.4ip . Moreover, suppose that 
(|6T24)1 - ((6T26|) hold on C x {n } x (t , r ~ T 0l t , r + T ). Then ([6^24]) - ([BT26]) /iold on 

C X Z X (f 0jJ . — T , ^0,r + 2o); 

F p j r (z,n,t r ) = 2(F p (z,n,t r )G r+1 (z,n,t r ) 

- G p+ i(z,n,t r )F r (z,n,t r )), (6.50) 
G p+Mr (z,n,i r ) = 4a(n,t r ) 2 (F p (z,n,t r )F+(z,n,t r ) 

-F+(z,n,t r )F r (z,n,t r )), (6.51) 
-R2 P +2(z) = G p+ i(z,n,t r ) 2 - 4a(n,t r ) 2 Fp(z,n,t r )F+(z,n,t r ), (6.52) 
(z, u,i r )eCxZx (t 0;T . — T , t ,r + To)- 

Moreover, 

4> tr (P,n,t r ) = -2a(n,t r )(F r (z,n,t r )<p(P,n,t r ) 2 + F+(z,n,t r )) 

+ 2{z - b + (n, t r ))F+{z, n, t r )4>{P, n, t r ) (6.53) 
+ (G> +1 (z,n,i r ) - G r+ i(z,n, t r ))(f>(P,n, t r ), 
at r (n, t r ) — — a(n, t r )(2(z — b + (n, t r ))F^~(z, n, t r ) 

+ G+ +1 (z,n,t r ) + G r+1 (z,n,t r )), (6.54) 
bt r (n, t r ) — 2((z — b(n, t r )) 2 F r (z, n, t r ) + (z — 6(n, t r ))G r +i(z, n, t r ) 

+ a 2 (n, t r )F^(z 7 n,t r ) — (a~(n, t r )) 2 TT(z, n,t r )), (6.55) 
(z,n,t r ) G C x Z x (i , r - T ,t .r +T ). 

Troo/. By LemmalEJwe have ([5T22]) . (|Q3j) . (|B37p . (j5^9|) - (j53T|l . and (j6T5V(l6T8l) 

for (n, i r ) G Z x (io. r — To, io. r - + T)) at our disposal. Differentiating (|6.52p at n = no 
with respect to t ri inserting (|6.50| and (|6.51|) at n = n , yields 

2T+a tr + aF„| tr = 2a(G p+ ii 7 ' J + - T+G r+ i) 

~ ~ ~ ~ ~ (6.56) 

= 2F+a( - 2{z - b+)F+ - G+ +1 - G r+1 ) + 2a(F+G+ +1 - G+ +1 F+) 

at n = no- By inspection, 

T+(z)G+ +1 (z) - G+ +1 {z)F+{z) = OOO. (6-57) 

|^| — *0O 

This can be shown directly using formulas such as (|2~2lI| - (l2T2"6"]) . (pT2"j) . (jO]) . (|63)l . 
and (|6.6p . It also follows from (|5.43p and the fact that F p is a monic polynomial of 
degree p. Thus one concludes that 

2F+a tr = 2F+a( - 2{z - b+)F+ - G+ +1 - G r+1 ) (6.58) 

at n = no, and upon cancelling 2T+ that (|6.54p holds at n = no. This and (|6.56|) 
then also proves that (|6.50p holds at n = no + 1- 

Next, differentiating 2aF p (f> = y — G p +i at n = no with respect to t r inserting 
(|6~50"1) . (j63T|) . and (|6~54]) at n = n , and using ([533]) to replace 2<zF+ by -(y + 
G p+ i)0 and (|5.22p to replace (G p+ i — y) by —2aF p (p, yields (|6.53p at n — n upon 
cancelling the factor 2aT p . 
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Differentiating (|6.46|) with respect to t r (fixing n = no), inserting (|6.46p (to 
replace Gp +1 ), (|6.51|) at n = no, and (|6.50|) at n = no + 1 yields 

= -2F+(b+ - 2(z - b+) 2 F+ + 2a 2 F r - 2(z ~ b+)G+ +1 ) 

+ 4(z - b+) 2 F+F+ + 4(z - b+)G p+1 F+ + 4(a) 2 F p F+ + G+ +Mr 
= -2F+(b+ - 2(z - b+) 2 F+ - 2(z - b+)G+ +1 + 2a 2 F r - 2(a+) 2 F++) 
- 4(a+) 2 F+F++ + A{z - b+) 2 F+F+ + 4(z - b+)G p+1 F+ 
+ 4a 2 F p F+ + G; +htr 
= -2F+(b+ - 2(z - b+) 2 F+ - 2(z - b+)G+ +l + 2a 2 F r - 2(a + ) 2 F++) (6.59) 
+ G+ +Mr - 4(a + ) 2 F+F+ + + (4a 2 F p + 4(z - b+) 2 F+ + 4(z - b + )G p+1 )F+ 
at n — no- Combining (|6.46|) and (16.47)) at n — hq one computes 

4(a+) 2 F++ = 4a 2 F p + 4(z - b+) 2 F+ + 4(z - b+)G p+1 (6.60) 
at n = no- Insertion of (|6. 60|) into (|6.59p then yields 

= -2F.+ (b+ - 2(z - b+) 2 F+ - 2{z - b+)G+ +1 + 2a 2 F r - 2(a+) 2 F++) 

+ G+ +htr 4(a+) 2 F+F++ + 4(a+) 2 F++F+ (6.61) 
at n — no- In close analogy to (|6.57|) one observes that 

F+(z)F+ + (z)-F++{z)F+{z) = 0(z p ^ 1 ) for peN. (6.62) 

\z\— »00 

Thus, since F^ is a monic polynomial of degree p, (|6.61[) proves that 

6+ - 2(z - b+) 2 F+ - 2{z - b+)G+ +1 + 2a 2 F r - 2(a+) 2 F++ = (6.63) 

at n — hq, upon cancelling . Thus, (|6.55|) holds at n = no + 1. Simultaneously, 
this proves (|6.5ip at n = no + 1. 

Iterating the arguments just presented (and performing the analogous consider- 
ations for ii < no) then extends these results to all lattice points n € Z and hence 
proves (|630] ) -(|535 ]) for (z, n, t r ) e C x Z x (t . r - T , t . r + T ). □ 

We summarize Lemmas I6J2J and I6JH next . 



Theorem 6.4. Assume Hypothesis 16.11 and condition (|6.41j) . Moreover, suppose 
that 

fj = fj{no,t r ), j = l,...,p, 
9j = 9j(n ,t r ), j = l,...,p-l, 

(6.64) 

9p + fp+l = 9p{ n 0, tr) + f P +i(t r ) 
for all t r e (t , r - T ,to, r + T ), 

satisfies the autonomous first-order system of ordinary differential equations (|6.27|) 
(for fixed n = no), 

fj.tr = Fj(f\, . . . , f p ,gi, . . . , g p -i, g p + f P +i), j = 1, . . . ,p, 

9j,t r =0j(fi,---,fp,9l,---,9p-i,9 P + fp+l), j = l,...,p-l, (6.65) 

(.9p + fp+i)t r = Gp{fi, ■ ■ ■ , f P , 9i, ■ ■ ■ , 9p-i,9p + fp+i) 
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with initial condition 

fj(n ,t . r ), j = l,...,p, 

9j{no,t , r ), j = l,...,p-l, (6.66) 
g p (n ,t 0r ) + / p +i(to,r)- 

Then F p and G p +i as constructed in (|6. 2(1 — f)6 .44() on C x Z x (t 0l , — To,to,r + To) 
satisfy the zero-curvature equations (|5.8[) . (15. 9ft . and (|5.45p . 

t/* r + UVr+i - Klj. X 17 = 0, (6.67) 

UV P+1 -V+ +1 U = 0, (6.68) 

V p+ i, tr -[V r+1 ,V p+1 }=0 (6.69) 
on C x Z x (io, r — Tq, to >r + ?b), 

with U, Vp+i, and V^+i given by (|5.10[) . In particular, a, 6 satisfy the algebro- 
geometric initial value problem (|5.2|) . (|5.3[) 



V b *r + 9p+i(a, b) - g +1 (a, 6) / 



K + 9 P +i (a,b)- g p+1 (a, b) J ( 6 70 ) 

M)U , = (« (0 ^ (0) )> 

on Z X (i ,r - ?0, *0,r + T ), 

and are given by 

tfnAO (^(fl 2p+2 ( Z )l/2)/^ 



a(n,t r ) 2 = i ^ 



lz=A»fc(™,*7-) 



2 ^ (Pfc(»,*r)-1)! 
q(n,t r ) 

x [] *0*fc(n > *r)-A**'(n,tr))~ , '* (n,tr) (6-72) 

k' = l, k'^k 

+ l(b^\n,t r )-b{n,t r f), 

^ 2p+l q(n,t r ) 

b{n,t r ) = - ^2 E m - ^2 Pk(n,t r )p k (n,t r ), (6.73) 

m=0 fc=l 

(z,n,t r ) e Z x (t 0)7 . - T ,t , r + T ). 
Moreover, Lemmas 13.2113.41 and 15.2115.51 apply. 



As in the stationary case, the theta function representations of a and b in the 
time-dependent context can be derived in complete analogy to the self-adjoint case. 
Since the final results are formally the same as in the self-adjoint case we again 
just refer, for instance, to [Sj, [7j, [S], [TO], [HI Sect. 1.4], [TH] (cf. also the appendix 
written in [5]), ES] , [301 Appendix, Sect. 9], [31 Sect. 13.2], [3J Sect. 4.6], [Ml Ch. 
28]. 

As in Lemma 14.21 we now show that also in the time-dependent case, most initial 
divisors are nice in the sense that the corresponding divisor trajectory stays away 
from infinity for all (n, t r ) S 1 x R. 
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Lemma 6.5. The set Mi of initial divisors D^ noto r \ for which I^/i( ni t r ), defined 
via (|5.53[) . is nonspecial and finite for all (n,t r ) eZxR, forms a dense set of full 
measure in the set Sym p (/C p ) of nonnegative divisors of degree p. 



Proof. Let A^ s ing be as introduced in the proof of Lemma I4~2l Then 



U (a Pa (M sing ) + t r U r 



U Upo ) + a Po (Sym^ 1 (K p )) + i r E/f ) (6.74) 



is of measure zero as well, since it is the image of 1 x Sym p_1 (/Cp) which misses 
one real dimension in comparison to the 2p real dimensions of J{JC P ). But then 

(J (^p (M)+nA Px _(P oo+ ) + t r fj^ ) ) (6.75) 

is also of measure zero. Applying a^ 1 to the complement of the set in (|6.75|) then 
yields a set .Mi of full measure in Sym p (/C p ). In particular, Mi is necessarily dense 
in Sym p (/C p ). □ 



Theorem 6.6. Let t 0r ) G Mi be an initial divisor as in Lemma 16.51 Then 

the sequences a,b constructed from //(no,to,r) «s described in Theorem 16.41 satisfy 
Hvvothesis \5 . 11 Ln particular, the solution a,b of the algebro- geometric initial value 
problem (|6T70|l . (loTTl]) is global in (n,t r ) eZxl. 

Proof. Starting with "Dft( no ,t r ) e •Mil the procedure outlined in this section and 
summarized in Theorem 16.41 leads to T>^ njt \ for all (n, t r ) eZx {to, r — T , t 0j ,+To) 
such that (|5.53p holds. But if a, b should blow up, then T>p / / ritt \ must hit Poo + which 
is impossible by our choice of initial condition. □ 

Note, however, that in general (i.e., unless one is, e.g., in the special periodic or 
self-adjoint case), "Dp.( n ,t r ) wu l S e t arbitrarily close to Poo + since straight motions 
on the torus are generically dense (see e.g. [H Sect. 51] or [TTJ Sects. 1.4, 1.5]) and 
hence no uniform bound on the sequences a(n,t r ),b(n,t r ) exists as (n,t r ) vary in 
ZxK. In particular, these complex-valued algebro-geometric solutions of the Toda 
hierarchy initial value problem, in general, will not be quasi-periodic (cf. the usual 
definition of quasi-periodic functions, e.g., in [311 p. 31]) with respect to n or t r . 

Appendix A. Hyperelliptic Curves of the Toda-Type 

We provide a brief summary of some of the fundamental notations needed from 
the theory of hyperelliptic Riemann surfaces. More details can be found in some of 
the standard textbooks [TT] and [57j as well as in monographs and surveys dedicated 
to integrable systems such as Ch. 2], [5], [H App. A, B], [31 App. A]. 

Fix p £ N. We intend to describe the hyperelliptic Riemann surface /C p of genus 
p of the Toda- type curve (|2.43[) . associated with the polynomial 

F P (z,y) = y 2 - R2 P +2(z) = o, 

2p+l (A 1) 

R2 P+ 2(Z) - J] ( Z £ ™)> {^}m=0 C C. 
m—0 



THE TODA HIERARCHY INITIAL VALUE PROBLEM 



3!) 



To simplify the discussion we will assume that the affine part of IC P is nonsingular, 
that is, we assume that 

E m ^ E m > for m^m', m, m! = 0, . . . , 2p + 1 (A. 2) 

throughout this appendix. Next we introduce an appropriate set of (nonintersect- 
ing) cuts Cj joining E m ^ and E m iij\, j = 1, . . . ,p + 1, and denote 

P+i 

C=\JC 3 , C j nC k = 0, j^k. (A.3) 
J=l 

Define the cut plane 

II = C\C, (A.4) 
and introduce the holomorphic function 

/2p+l s 1/2 

i? 2p+2 (-) 1/2 : n^C, ( ]J(z-E m )) (A.5) 

^ m=0 ' 

on II with an appropriate choice of the square root branch in (|A.5|) . Next we define 
the set 

M p = {{z,aR 2p+2 {z) 1/2 ) | zeC, a E {1, -1}} U {P^, P^} (A.6) 

by extending i?2p+2( • ) 1 ^ 2 to C. The hyperelliptic curve K, p is then the set M. p with 
its natural complex structure obtained upon gluing the two sheets of M. p crosswise 
along the cuts. Moreover, we introduce the set of branch points 

B{K, p )^{{E m ^tl. (A.7) 

Points P E IC P \ {Poo± } are denoted by 

P=(z,aR 2p+2 (z) 1 / 2 ) = (z 1 y), (A.8) 

where y(-) denotes the meromorphic function on K p satisfying J- p {z,y) = y 2 — 
P2 P +2(z) = and 

y{P) f = Q t(i - \ ( £ E ™) C + °(C 2 )) r p " x as P - P 0O± , c = 1A. (A.9) 

In addition, we introduce the holomorphic sheet exchange map (involution) 
*:K p ^K Pl P=(z,y)»P* = (z,-y), P 00± ^P^ ± =P 00 ^ (A.10) 
and the two meromorphic projection maps 

7r: IC p — > C U {oo}, P = (z, j/) i— > z, Poo ± !— > oo (A. 11) 

and 

j/: £ P -^CUH P=(z,y)^y, P QO± ^ oo. (A.12) 

Thus the map ir has a pole of order 1 at Poc ± and y has a pole of order p + 1 at 
Poo ± . Moreover, 

= tt(P), y(P*) = -y(P), PE/C P . (A.13) 

As a result, K, p is a two-sheeted branched covering of the Riemann sphere CP 1 
(= C U {oo}) branched at the 2p + 4 points {(E m , O)} 2 ^^ , P OQ± . K p is compact 
since tt is open and CP 1 is compact. Therefore, the compact hyperelliptic Riemann 
surface resulting in this manner has topological genus p. 



40 



F. GESZTESY, H. HOLDEN, AND G. TESCHL 



Next we introduce the upper and lower sheets IT± by 

n± = {(z, ±R 2p+2 (z) 1/2 ) e M p | z e n} (A.u) 

and the associated charts 

C±:H ± ^n, P^z. (A.15) 

Let {dj,bj}^ =1 be a homology basis for JC P with intersection matrix of the cycles 
satisfying 

a j obk=5j ! k, a,joak = 0, b ] ob k =0, j,k = l,...,p. (A. 16) 

Associated with the homology basis {a,j, bj}^ =1 we also recall the canonical dis- 
section of IC P along its cycles yielding the simply connected interior IC P of the 
fundamental polygon dK. p given by 

dlC p = aibia^b^ a^bia^, b^ ■ ■ ■ a^b^ 1 . (A. 17) 

Let M(IC n ) and A / ( 1 (/C n ) denote the set of meromorphic functions (0- forms) 
and meromorphic differentials (1-forms) on K, n . The residue of a meromorphic 
differential v G A4 1 (/C„) at a point Q G JC n is defined by 

res Q (,) = ^/^, (A.18) 

where 7q is a counterclockwise oriented smooth simple closed contour encircling Q 
but no other pole of v. Holomorphic differentials are also called Abelian differentials 
of the first kind. Abelian differentials of the second kind a/ 2 ) G A4 1 (/C„) are 
characterized by the property that all their residues vanish. They will usually be 
normalized by demanding that all their a-periods vanish, that is, 



/ w (2) =0, j = l,..., 

J a 4 



n. (A.19) 



(2) p 

If Wp n is a differential of the second kind on K, n whose only pole is Pi £ IC n with 
principal part £ _n_2 d(, n G No near Pi and uij = ( J2m=o dj,m{Pi)( m ) d( near Pi, 
then 

ii^ = %TT- "-M.... (A.20, 

Using local charts one infers that dz/y is a holomorphic differential on K p with 
zeros of order p — 1 at Poo ± and hence 

zi~ x dz , . 

Vj = ——, j = l,..., p, (A.21) 

form a basis for the space of holomorphic differentials on /C p . Introducing the 
invertible matrix C in C p 

C ={ C i*) jt k=x,...j C iM = I '/.- (A.22) 
c(fc) = ( Cl (fc),...,c p (fc)), Cj (fc) = (C" 1 ) j,fc = l,...,p, (A.23) 

the normalized differentials ojj for j = 1, . . . ,p, 

p , 

Uj = Cj(t)Vi, / = 6 jtk , j,k = l,...,p, (A.24) 

J ah 
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form a canonical basis for the space of holomorphic differentials on IC P . 
In the chart {Up 00± ,(p 00± ) induced by 1/tt near P oc± one infers, 



p 



"-^■■■•■^-^tt^Ftfw <a ' 25) 



= ±^(p) + cQc(p) £ £ m + c(p-l)) +0(C 2 ))^C as P-P^, 

c = v*. 

The matrix r = (t 7 ^)^_ 1 in C pxp of ^-periods defined by 

Tj,e= [we, j,£ = l,...,p, (A.26) 
Jb 3 

satisfies 

Im(r) > and T jtt = n tj , j, £=l,...,p. (A.27) 
Associated with the matrix r one introduces the period lattice 

L p = {zeC p \z = m + nr,m,neZ p } (A.28) 
and the Riemann theta function associated with tC n and the given homology basis 

{dj, &j}j=l,...,n; 

0(z) = ^2 cxp (27ri(n,z) +7ri(n,nr)), zeC", (A.29) 

where = ot t = Y^j=i^j v j denotes the scalar product in C". It has the 

fundamental properties 

. . . . .,z n ) = 6(z), (A.30) 

0(z + m + nr) = exp ( - 2m{n,z) - m(n,nr))6(z), m,neZ n . (A. 31) 

Next, fixing a base point Qo € K, p \ {Poo±}, one denotes by J(fC p ) = C p /L p the 
Jacobi variety of JC P , and defines the Abel map Aq by 

Aq ( /- £n ^ J{lCp), A Qo (P)=( / uji,..., u> p ) (mod L p ), P e K, p . 

\jQo JQo / 

(A.32) 

Similarly, one introduces 

a Qo : Div(/Cp) - J(/C p ), V ~ a Qo (2?) = £ P(P)A Qo (P), (A.33) 

where Div(/C p ) denotes the set of divisors on JC p . Here a map V: JC p — > Z is called 
a divisor on /C p if 2?(P) ^ for only finitely many P G JC p . (In the main body of 
this paper we will choose Qq to be one of the branch points, i.e., Qq € B{K, P ), and 
for simplicity we will always choose the same path of integration from Q to P in 
all Abelian integrals.) 

In connection with divisors on IC P we will employ the following (additive) nota- 
tion, 

V Qo q = V Qo +Vq, Vq = V Qi +---+ V Qm , (A.34) 
Q = {Qi, Q m } G Sym m K p , Qa G JC p , me N, 
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where for any Q £ K, p , 

^ Kp -n„, p^ o( p) = {; 'Zp^m < A ' 35 » 

and Sym m K, p denotes the mth symmetric product of K, p . In particular, Sym m K, p 
can be identified with the set of nonnegative divisors < T> £ Div(/C p ) of degree m £ 

N. A divisor V Q = V Ql H h 2? Qm will be called finite if Q k £ /C p \ {Poo+ , -Poo_ }, 

fc = 1, . . . , m. 

For / £ M{K, P ) \ {0}, a; £ Al 1 (/C p ) \ {0} the divisors of / and oj are denoted 
by (/) and (u), respectively. Two divisors V, £ £ Div(/C p ) are called equivalent, 
denoted by V ~ £, if and only if V - £ = (/) for some / £ M(JC p ) \ {0}. The 
divisor class [V] of V is then given by [V] = {£ £ Div(/C p ) | £ ~ 2?}. We recall that 

deg((/))=0, deg((o;))=2(p-l), feM(JC p )\{0}, u; E M 1 (IC P ) \ {0} , 

(A.36) 

where the degree deg(X>) of V is given by deg(2?) = J2pek p ^ ^ s customary 

to call (/) (respectively, (ui)) a principal (respectively, canonical) divisor. 
Introducing the complex linear spaces 

C(V) = {/ £ M(K p ) | / = or (/) > V}, r(V) = dim c C(V), (A.37) 

C\V) = {loE M\IC P ) \lo = 0ot(lo)> V}, i(V) = dim c C X (V) (A.38) 

(with i(T>) the index of specialty of V), one infers that deg(X>), r(T>), and i(D) only 
depend on the divisor class [V] of V. Moreover, we recall the following fundamental 
facts. 

Theorem A.l. Let V £ Div(/C p ), w £ .M 1 ^) \ {0}. T/ien, 

i(£>) = r(£> - (w)), v G N . (A.39) 
TTie Riemann-Roch theorem reads 

r(-V)=deg(V)+i(V)-p+l, n £ N . (A.40) 
-By Abel's theorem, V £ Div(£ p ), p £ N, is principal if and only if 

deg(V) = and a Qo (2?) = 0. (A.41) 

Finally, assume p £ N. TTien ag o : Div(/C p ) — > J(JC P ) is surjective (Jacobi's 
inversion theorem). 

Theorem A.2. Let V Q £ Sym p K p , Q = {Q 1 ,..., Q p }. Then, 

1 < %{Vq) = s (A.42) 

if and only if there are s pairs of the type {P, P*} C {Q\, . . . , Q p } (this includes, 
of course, branch points for which P = P*). Obviously, one has s < p/2. 

Next, we denote by Eq o = (Eq x , . . . , Sq 0jP ) the vector of Riemann constants, 

1 P f f P 

2q„, 3 = 2(1 + T jd ) - / w/(P) / Uj, j = l,..., P- (A.43) 

l— l J at •> Qo 

Theorem A. 3. Let Q = {Qi, ■ ■ ■ , Q p } £ Sym p IC P and assume T>q to be nonspecial, 
that is, i(T>Q) = 0. Then, 

0(Z Qo - A Qo (P) + a Qo (Dg)) = if and only if P £ {Q u . . . ,Q p }. (A.44) 
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Appendix B. Some Interpolation Formulas 



In this appendix we recall a useful interpolation formula which goes beyond the 
standard Lagrange interpolation formula for polynomials in the sense that the zeros 
of the interpolating polynomial need not be distinct. 

Lemma B.l. Let p£N and S p -i be a polynomial of degree p — 1. In addition, let 
F p be a monic polynomial of degree p of the form 



F P (z) =t[(z »k) Pk , Pj G N, fij e C, j = 1, 



k=l 



,q, ^2pk=p- 

fc=i 



Then. 



t\{p k -i-l)\ 



fe=l 1=0 



Pk - 



fc'=i, fc'^fc / 



In particular, S p -\ is uniquely determined by prescribing the p values 

ai i/ie given points fx%. . . . , fj, q . 

Conversely, prescribing the p complex numbers 



(o) (i) 



k ' u fc i • • • J u fc 



„(Pfc-l) 



fc = 1,. 



,q, 



there exists a unique polynomial T p _i of degree p — 1, 



q Pk~ 1 

T p ^(z)=F p (z)J2J2 



a 



k=l i=Q 



e\( Pk -£-iy. 



such that 



T p -i(fJ,k) = a£\T;_M = af, . . .,T™- lJ (»k) = ^ 



Pk-t- 



i i^-cr 1 n (c-w* 



-Pfc' 



(i) 



(B.l) 
(B.2) 

zeC. 

(B.3) 

(B.4) 
(B.5) 

zeC, 



fc = 1, . . . ,q. 



(B.6) 

Proof. Our starting point for proving (|B.2[) is the following formula derived, for 
instance, by Markushevich [H Part 2, Sect. 2.11, p. 68], 



Sp-i(z) 



1 

2~7ri 



dCS p ^(0F p (C)-F p (z) 



F P (C) 



Z £ 



(B.7) 



'r x pv 

where T is a simple, smooth, counterclockwise oriented curve encircling . . . , fi q 
strictly in its interior. Since the integrand in (|B.7|) is analytic at the point C, — z, 
we may, without loss of generality, assume that Y does not encircle z. With this 
assumption one obtains 

i r dCSp-iiO 

2tt« J r C - z 







(B. 
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and hence deforming T into sufficiently small counterclockwise oriented circles Tk 
with center at /ik, k = 1, . . . ,q, such that no \xy , k' ^ k, is encircled by Tk, one 
obtains 



Sp-i(z) 



F p (z) I dC5 p _i(0 



2m J r F p (()((-z) 
F P (z) I dC-gp-i(C) 



F P (()(C-z) 




rfC(C-Mfc) 




(C-*) n^C- 




dC (C - 


"Pfc 


(c-*)ii£ =1 (c- 


fl k ' ) Pk ' 






dC(C-^) £ 


-Pk 



2m l\ 

k=i i=o 

2ni 2-^2-^ £\ 

k=l 1=0 

Fpjz) y^y^ 1 Sp-liP-k) I <n_ K -- in,- I ' " 



2m tx to el ^ k - *) ift =1 (c - 

where we used 

' dC (C - Mfc)'- p */(C) = if £ > Pk, * G N, (B.IO) 

r fc 

for any function / analytic in a neighborhood of the disk Dk with boundary Tk, 
k = 1 , . . . , q, to arrive at the last line of (|B.9j) . An application of Cauchy's formula 
for derivatives of analytic functions to (|B.9[) then yields 

k=l t=o 

x *d£ k dC (C - w)^-'-^ 1 (C - z) Ul= hk ^k(C ^ 

q Pfc-i q(0 
p- 



= WEE ' s;j|,/ ' i ' i 



d Pfe " 



, zeC, 

(B.ll) 



and hence (|B.2[) . Conversely, a linear algebraic argument shows that any polynomial 
T p _i of degree p — 1 is uniquely determined by data of the type 

Tp^i^^M,...,^-^^), k = l,...,q. (B.12) 
Uniqueness of the representation (|B.2p then proves (|B.5|) . □ 
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We briefly mention two special cases of (|B.2j) . First, assume the generic case 
where all zeros of F p are distinct, that is, 

q = p, Pk = 1, fJ>k 7^ Mfc' for k ^ k', k,k' = 1, ...,p. (B.13) 
In this case (|B.2I) reduces to the classical Lagrange interpolation formula 

^'''° w £ 'woMil"(.- W ) ' ,B - 14) 

Second, we consider the other extreme case where all zeros of F p coincide, that is, 
q = l, Pl =p, F p (z) = (z- Aii) p , zeC. (B.15) 
In this case (|B.2[) reduces of course to the Taylor expansion of S p -i around z = /ii 

p-i S W / \ 

g p -i(^=E p\ (z-ViY, z^C. (B.16) 

Appendix C. Asymptotic Spectral Parameter Expansions and 

Nonlinear Recursion Relations 

In this appendix we discuss asymptotic spectral parameter expansions for F p /y 
and Gp+i/y as well as nonlinear recursion relations for the corresponding homoge- 
neous coefficients fg and gi and analogous quantities fundamental to the polynomial 
recursion formalism for the Toda hierarchy. 

We start by recalling the following elementary results (which are consequences 
of the binomial expansion) . Let 

{£ ; m } m =o,...,2 P +i C C for some p E N (C.l) 

and r\ G C such that \r\\ < min{|i?o| _1 , • ■ ■ , |i?2p+i 

Then 

• 2p+l . -1/2 oo 

n(i--M) =£4(£>A (c.2) 

m=0 ' fc=0 

where 

c Q (E) = 1, 

o k {E)= t Mj0t^ E °--- E ^> (C - 3) 

JoH hj2p+i=fc 

The first few coefficients explicitly read 

1 2p+l 

co(s) = i, oiCD - -^ £ m , 

m=0 



£a(^) = 2 £ « 3 + jE^' etc - (°- 4 ) 

Similarly, 



4 

mi, 7712—0 771 — 

7711 <7Tl 2 



2p+l \ 1/2 OO 

m=0 ' k=0 
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where 
c (E) = 1, 

(2j y....{2j 2p+1 y.Et ■■■e%+± 



, i _ n 2 (Jo!) •■ ■(j 2 p+i!) 2 (2jo - 1) •■ -(2^+1 - 1) 



ioH K?2p+l=fc 

The first few coefficients are given explicitly by 

1 2p+l 

Co® = 1, C!(E) =-^J2 Em > 



2 

m=0 



2p+l 2p+l 

c 2 (E) = J E E mi E m2 --Y,E 2 ml etc. (C.7) 



4 

mi ,7712—0 
777,1 <777 2 



Theorem C.l. Assume (|2.1[) . s-Tl p (a, 6) = 0, and suppose P — (z,y) € K p \ 
{P 00+ ,P 00 _}. T/ien Pp/y and Gp+i/y have the following convergent expansions as 
P — > P 

y *!=0 y £=-1 

where £ = 1/z is i/ie /oca? coordinate near P oc± and and gi are the homogeneous 
versions of the coefficients ft and ge introduced in (|2.8[) . In particular, fi and gi 
can be computed from the nonlinear recursion relations 

/o = l, fi = -b, f 2 =a 2 + (a-) 2 + b 2 , 
^ t+i e+i 

fe+2 — ~-z E ft+2-kfk - ft+l-kfk 
k=l fe=0 

e-i 

+ E ( " 26 3 A-i-fe/fe + 2a 2 bf+^_J k + 2(a~) 2 bf e - 1 _ k f-) 



k=0 
1-2 



E {^ 2 ft-2-Jk + {a-) 2 b 2 U-2-kf k - + a 2 (a-) 2 //_ 2 - fc A- 



fe=0 

-^ 4 //_2-^-^(0 4 /,-2- fc /^) ) ^N. (C9) 



and 



.9-1 = -1, .90 = 0, gi = -2a 2 , 

1 1 I 1 

9e+i = 5 E ( 6 + b+ )k-i-kgk + r 22,k-k9k (CIO) 



fc=-l fc=0 



1 

o E ( bb+ 9t-2-k9k - a 2 (.9£_ 2 -fc +9i-2-k)(9k + §t ))> ^N. 



fc=-i 
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Moreover, one infers for the E m - dependent summation constants ct, I = 0, . . . ,p+l, 
in F p and G p +i that 

c e = c e (E), t = 0,...,p+l (C.ll) 

I 

fe = J2 c i-k(K)f k> e = 0,...,p, (C.12) 

fc=0 



9i + fp+i$ P ,e = ^2ce-k(R )g k - c e+1 (E ), i = Q,...,p, (C.13) 

tAp 

= X>_ fe (£)/*, eeN 0l (C.14) 



fc=0 
lAp 



91 



^de-k(E)(9k + fp+iS P ,k)-£e+i(E), £eN . (C.15) 

fe=0 

id Gp+i b 

z p+ 1 near infinity) one obtains 

Fp[z) (^cm*- k ){j:hz- l - 1 )=Y.fc- i -\ (c.i6) 



Proof. Dividing F p and G p +i by R^ p 2 +2 (temporarily fixing the branch of i?2p+2 as 



R2p+2(z) 1 / 2 |*Hoc 

G P+l( z ) 
R2p+2{z) 1 / 2 khoo 



p+1 \ oo 



E4©^ E^' =^ E ( c - 17 ) 



for some coefficients and ^ to be determined next. Here we have temporarily 
introduced the notation 

p p+i 
G p+1 (z) = -z? +1 + £ 9p-d + f P+1 = 9v-t*. (C18) 

Dividing (|2.37[) and (|2.39p by i?2p+2 and inserting expansions (|C.16[) and (|C17() 
into the resulting equations then yield the nonlinear recursion relations (|C.9[) and 
(|C.10P (with ft and gi replaced by fi and gt, respectively). More precisely, one 
first obtains |/o| = \g~i\ = 1 and upon choosing the signs of fo and <?_i such that 
fo = fo = 1 and g_ x = -1 one obtains (|C.9|) and (|C.10j) . Next, dividing (|23T|) and 
(|2.32[) by i?2p+2' inserting the expansions (|C. 1 6|) and (|C.17|) . and comparing powers 
of z~ e as z — > oo, one infers that /g and satisfy the linear recursion relations 
(|2.4p - (|2.6p . Hence one concludes that 

// = //, 91= 9l, ^eN (C.19) 

for certain values of the summation constants q. To show that ft — fi, gg = gt, 
and hence all q, £ £ N, vanish, we recall the notion of degree as used in the proof 
of Lemma 15.41 which serves as an efficient tool to distinguish between homogeneous 
and nonhomogeneous quantities. To this end we employ the notation 

f(r) =S (r) f) / = {/(„)} neZ C C, 5« = (5^ r _'r '"n r G Z ' ( C ' 2 °) 

\(S ) r , r < 0, 



m A 71 = min{m, n}. 
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and introduce 

deg(a (r) ) = deg(6 (r) ) = 1, r e Z. (C.21) 

This results in 

deg{f t )=l, deg(g i )=£ + l, leN. (C.22) 
using induction in the linear recursion relations (|2.4p - (|2.6j) . Similarly, the nonlinear 
recursion relations (|C9|) and (IC.10|) yield inductively that 

deg(fi)=£, deg(g e )=£ + l, feN. (C.23) 
Hence one concludes that 

h = h, 9t = 9h ^N„. (C.24) 

A comparison of coefficients in (|C.f 6|) proves (|C.14|) . Similarly, we use (|C.17[) to 
establish ()C15j) . Next, multiplying (|C.2j) and (|C.5j) . a comparison of coefficients of 
rj k yields 

k 

Ck-i(E)c £ (E) = 4.o, k g N . (C.25) 

e=o 

Thus, one computes 

l I m It 

^ cz- m (E_)f m = ^ ^ ce- m (E)c m - k (E)f k = YYci^ p (E_)c p ^k(E_)f k 

m=0 m=0 k=0 k=0p=k 

i , i-k v 

= ( S «_ fc _ m (£)c m (E) J/ fc = f e , £ = 0,... ,p, (C.26) 

fc=0 ^m=0 ' 

applying ([C?25)) . Hence one obtains ([CT2]) and thus ([CTT]) (cf. ((2U)). The corre- 
sponding proof of (|C.13[) is similar to that of fi. □ 
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